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The  hugoniot  equations  of  state  have  been  determined  experi¬ 
mentally  for  the  materials  over  the  following  pressure  ranges: 
Fans teel-77  (0.3  to  5.0  Mb),  OFHC  copper  (1.0  to  4.5  Mb),  NTS 
Climax  Stock  granodiorite  (0.6  to  1.9  Mb),  NTS  Rainier  Mesa 
tuff  -  vacuum  dry  -  (0.2  to  1.3  Mb) ,  NTS  Rainier  Mesa  tuff  - 
>  95%  water  saturated  -  (0.08  to  1.33  Mb),  V44Q  grouting  mate¬ 
rial  (0.11  to  1.15  Mb),  schistose  gneiss  (0.^7  to  2.5  Mb), 

NTS  Ferris  Wheel  dolomite  (0.63  to  1.74  Mb),  NTS  Banded  Moun¬ 
tain  limestone  (0.15  to  1.22  Mb),  NTS  desert  alluvium  (0.15 
to  0.90  Mb).  Where  possible,  comparisons  are  made  to  previous 
hugoniot  measurements  by  other  workers . 
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FOREWORD 

This  report  presents  a  description  and  gives  the  results  of 
the  experimental  work  performed  under  contract  DA-49-146-XZ- 
429.  The  work  was  sponsored  by  the  Strategic  Structures 
Division  of  the  Defense  Atomic  Support  Agency,  Washington, 

D.  C.,  during  the  period  October  1965  to  December  1968.  High 
pressure  hugoniot  state  measurements  were  performed  on  five 
Nevada  Test  Site  (NTS)  rocks,  a  schistose  gneiss  and  a 
volcanic  tuff  impedance-matching  grout.  In  addition  several 
experiments  were  performed  or.  OFHC  copper  and  Fansteel-77, 
a  copper-nickel-tungsten  alloy.  These  two  metals  were 
employed  throughout  the  program  as  hugoniot  impedance  match 
standards.  Table  XIII  lists  the  materials  tested  and  the 
pressure  ranges  achieved. 

\ 

Section  I  of  the  report  describes  in  detail  the  experimental 
procedures  and  target  designs  employed^  for  the  hugoniot  deter¬ 
minations.  The  data  analysis  procedure  is  presented  in  Section 
II,  as  well  as  a  discussion  of  sources  of  experimental  error. 
Section  III  presents  the  experimental  data  and  descriptive 
interpretation  of  the  resulting  hugoniots.  No  attempt 
is  made  to  relate  the  results  to  geophysical  problems.  To 
ascertain  the  accuracy  and  validity  of  the  compiled  results, 
comparisons  are  made  with  the  published  data  of  other  re¬ 
searchers.  A  summary  of  the  results  is  presented  in  Section 
IV  and  general  conclusions  are  drawn. 
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INTRODUCTION 


The  purpose  of  the  program  described  herein  was  to  experi¬ 
mentally  obtain  high  pressure  shock  wave  data  to  aid  in  the 
interpretation  of  the  results  of  underground  nuclear  testing. 
Since  the  nuclear  experiments  are  performed  at  NTS,  the  inter¬ 
est  is  mainly  in  the  NTS  geological  environment.  The  experi¬ 
mental  program  undertaken  at  first  by  GM  Defense  Research  . 
Laboratories,  Santa  Barbara,  California,  and  then  the  Materials 
and  Structures  Laboratory  of  Manufacturing  Development,  General 
Motors  Technical  Center,  Warren,  Michigan,  have  been  to  provide 
data  on  several  selected  rocks  from  NTS. 

In  the  past  two  decades,  high  explosives  have  been  the  main 
energy  sources  for  experimental  systems  used  to  induce  high  pres¬ 
sure  shock  states  in  materials  of  interest.  By  these  techniques, 
pressures  of  about  1  Mb  may  be  achieved  in  the  denser  rocks .  For 
this  study  an  accelerated-reservoir  light-gas  gun  (ARLG  gun)  was 
employed  as  an  impactor  launching  instrument. ^  The  ARLG  gun 
is  capable  of  launching  high  density  metal  impactors  at  veloci¬ 
ties  of  up  to  ^  8  km/sec.  This  method  of  experimentation  offers 
significant  advantages  over  the  explosive  techniques  previously 
used.  Shock  induced  pressures  of  over  2  Mb  in  many  rocks  are 
obtainable  with  impactors  launched  by  the  ARLG  gun.  By  proper 
selection  of  impactor  material  and  launch  velocity,  the  induced 
shock  pressure  states  are  continuously  variable  over  the  range 
from  ^  0.1  to  2.0  Mb. 

The  experimental  program  described  here  may  be  considered  for 
convenience  to  consist  of  three  phases,  though  not  necessarily 
in  the  sequence  presented.  The  first  effort  was  to  establish 
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the  hugoniots  of  OFHC  copper  and  Fansteel-77.  The  second  phase 
involved  the  determination  of  the  hugoniot  equations  of  state 
for  NTS  Climax  Stock  granodiorite  ("granite") ,  Rainier  Mesa 
volcanic  tuff  (dry) ,  a  NTS  Ferris  Wheel  dolomite,  a  NTS  Banded 
Mountain  limestone,  and  a  schistose  gneiss,  not  from  NTS.  The 
third  phase  consisted  of  the  hugoniot  determination  of  materials 
that  presented  peculiar  experimental  difficulties,  including 
the  water-saturated  Rainier  Mesa  tuff,  the  grouting  material, 
V44d,  and  the  NTS  desert  alluvium. 
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SECTION  I 


HUGONIOT  STATE  MEASUREMENTS 


TECHNICAL  APPROACH 


Hugoniot  Theory 

Application  of  the  principles  of  conservation  of  mass, 

momentum,  and  energy  across  a  shock  front  have  led  to  the 

well-known  Rankine-Hugoniot  equations .  The  equations  may  be 

used  to  relate  the  discontinuous  changes  across  a  shock 

front  of  the  variables:  Pressure,  P,  density,  p,  specific 

volume,  V  =  — ,  particle  velocity,  u  ,  and  specific  internal 
P  P 

energy,  E.  Where  the  shock  wave  is  travelling  at  constant 

velocity,  U  illustrated  in  Figure  1.  The  particle  velocity 
s 

ahead  of  the  shock  is  zero,  and  subscript  1  designates 
conditions  behind  the  shock. 


p  U  =  p,  (U  -u  ) 
o  s  1  s  p 


P1  -  Po  -  po  Vp  (2) 

E1  -  Eo  =  p/2  (P1  +  V  <Vo  -  Vl>  (3) 

A  measurement  of  any  two  parameters,  e.g.  shock  wave  velocity 
and  the  particle  velocity,  associated  with  a  shock  wave  provides 
sufficient  information  to  determine  the  other  adiabatically 
achieved  state  variables  of  the  shocked  specimen.  implicit 
in  the  Rankine-Hugoniot  relations  is  the  assumption  that  the 


E. 
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material  behind  the  shock  front  has  reached  steady-state  or 
equilibrium  conditions. 


Material  At  Material  At 

Shock  Conditions:  Initial  Conditions . 


=  =C> 

Shock 

Front 

Velocity  -  Us 

Figure  1  Schematic  of  a  Shock  Wave  Propagating 
Into  a  Material 


Although  measurement  of  shock  velocity  is  relatively  straight¬ 
forward,  the  measurement  of  particle  velocity  is  more  difficult 

experimentally.  For  this  study  an  adaptation  of  the  impedance 

.  (2) 
matching  technique  developed  by  Walsh,  et.  al.  was  used  to 

calculate  the  particle  velocity  associated  with  a  measured  shock 

wave  velocity  in  the  specimen.  A  shock  wave  was  induced  in  a 

specimen  material  customarily  by  the  arrival  of  an  impactor  at 

its  front  surface.  The  shock  wave  then  passed  through  the 


specimen,  where  the  specimen  shock  wave  transit  time  was  recorded. 

The  shock  velocity  was  calculated  and  in  the  P-u^  plane  a  line 

of  slope  p  U  established  the  Rayleigh  line  passing  through 
o  s 

P  =  0,  Up  =  0,  as  shown  in  Figure  2. 
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Figure  2  Schematic  of  Impedance  Match  Method 
Using  Dissimilar  Materials 

The  hugoniot  of  the  impactor  in  the  form  P  =  P(u^)  previously 
determined  then  was  reflected  about  the  vertical  established 
by  ^  =  v  (the  impact  velocity) .  The  intersection  of  the 
Rayleigh  line  with  the  reflected  impactor  hugoniot  represents 
the  shock  state  of  the  specimen. 

An  advantage  of  substantial  importance  is  the  possibility  of 

impacting  the  specimen  with  a  gun  launched  impactor  of  the  same 

material.  With  this  technique ,  the  impedance  symmetry  is 

exact,  i.e.,  u  =  1/2  v.  The  shock  pressure  may  be  determined 
P 

directly  from  equation  (2) . 


Figure  3  illustrates  schematically  the  impedance  symmetry 
approach  wherein  like  materials  are  made  to  collide.  A 
series  of  experiments  on  OFHC  copper  and  Fansteel-77  were 
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performed  over  the  full  velocity  range  of  the  gun  for  each 
material  using  these  symmetrical  impact  conditions,  to 
determine  their  hugoniots.  The  accuracy  and  precision  with 
which  the  impact  velocity  may  be  measured  is  0.03%. 


Figure  3  Schematic  of  Symmetrical  Impedance 

Match  Method  t 

With  the  standard  hugoniots  established,  the  second  phase 
test  series  began  using  impactors  of  either  OFHC  copper  or 
Fansteel~77  and  the  dissimilar  impedance  approach  as  discussed 

above . 
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EXPERIMENTAL  TECHNIQUES 


During  the  course  of  the  work  on  the  program,  refinements  of 
design  and  instrumentation  have  been  made.  The  following  is 
a  description  of  the  experimental  system  as  it  is  currently  in 
use.  The  areas  of  significant  improvement  are  in  target  design 
and  in  the  instrumentation  for  the  measurement  of  impact  velo¬ 
city.  The  changes  will  be  noted  in  the  discussion  of  these 
topics . 


Launching  Techniques 

The  gun  used  to  accelerate  the  impactor  is  an  accelerated- 
reservoir  light-gas  gun  with  a  launch  tube  bore  diameter  of 
either  29  mm  or  64  mm.  This  type  of  gun  maintains  a  reasonably 
constant  pressure  on  the  base  of  the  projectile  during  the 
launch  allowing  the  acceleration  of  a  variety  of  unshocked 
impactor  materials.  Figure  4  shows  the  layout  of  the  range. 

The  gun  consists  of  the  following  major  components: 

1.  Powder  chamber 

2.  Pump  tube,  87  mm  inside  diameter  by  12  m  long 

3.  Accelerated-reservoir  high-pressure  coupling 

4.  Launch  tube,  either  29  mm  inside  diameter  by  8  m 
long  or  64  mm  inside  diameter  and  8  m  long 

5.  Instrumented  target  chamber  and  flight  range. 

When  the  gun  is  loaded  for  firing,  a  weighted  plastic-nosed 
piston  is  placed  in  the  breech  of  the  pump  tube.  Gun  powder 
is  then  loaded  in  to  the  powder  chamber  behind  the  piston 
and  the  pump  tube  is  filled  to  the  required  pressure  with 
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hydrogen.  After  ignition  of  the  gun  powder,  the  hydrogen  is 
compressed  into  the  coupling  by  the  accelerated  piston.  The 
saboted  impactor  is  in  turn  accelerated  by  the  release  of  the 
compressed  hydrogen  through  a  high  pressure  burst  diaphragm. 


I  V. 


Figure  4  Layout  of  ARLG  Gun  Range 

Prior  to  firing,  the  flight  range  and  instrument  chamber  are 
evacuated  and  then  flushed  with  helium  to  approximately 
10  3  torr  to  eliminate  any  spurious  effects  of  gas  build-up 
and  ionization  between  projectile  and  target.  Sealing  lips 
on  the  rear  of  the  plastic  sabot  are  pressed  tightly  against 
the  inside  of  the  launch  tube  bore  by  an  interference  fit 
and  the  high  pressure  gas,  effectively  eliminating  blow-by. 

Careful  attention  to  the  condition  of  the  launch  tube  is  nec¬ 
essary  for  successful  firing  in  the  high  velocity  ranges.  Bore 
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linearity  of  better  than  0.2  mm  over  the  full  8  m  length  of  the 
launch  tube  is  maintained.  The  internal  diameter  is  maintained 
to  within  0.01  mm.  Launch  tubes  are  cleaned  and  honed  after 
each  firing  and  are  removed  every  15  to  20  firings  for  thorough 
reconditioning . 

Figures  5  and  6  show  the  instrumentation  chamber  designed  for 
the  high  pressure  studies.  The  chamber  is  connected  to  the 
barrel  of  the  gun  only  through  an  0-ring  seal  to  allow  free 
axial  movement  of  the  launch  tube.  The  target  chamber  is 
shock-mounted  to  prevent  target  displacement  before  projectile 
impact.  To  facilitate  this,  several  stages  of  mechanical 
isolation  have  been  arranged  in  the  barrel,  I-beam  support 
structure  and  concrete  foundation. 


Figure  5  Target  Chamber  Set-up 
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STATION 

#1  STATION 


Figure  6  Schematic  of  Target  Chamber  (Top  View) 

The  target  chamber  is  a  steel  cylinder  of  61  cm  O.D.  and  1.5  m 
length.  Physical  access  and  instrument  ports  are  precisely 
machined  in  a  horizontal  plane  and  in  planes  45°  above  horizontal 

Two  sets  of  the  ports  are  employed  as  x-ray  windows  and  house 
the  secondary  spatial  references .  The  references  are  cross 
wires  with  the  separation  distance  between  windows  (nominally 
30.0  cm)  accurately  determined  by  an  optical  comparator. 

Operationally,  the  ports  are  closed  against  O-ring  vacuum  seals 
with  Plexiglas  or  magnesium  windows  for  optical  and  x-ray  instru¬ 
mentation  or  with  steel  cover  plates. 
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Instrumentation 

The  impactor  velocity  measuring  system  consists  of  a  two  laser 
triggering  system  and  a  short  duration  dual  flash  x-ray  as 
described  by  R.  Lingle  and  A.  H.  Jones. With  such  a  system, 
impactor  velocities  are  measured  accurately  to  'v  0.03%.  The 
triggering  system  consists  of  two  sets  of  neon-helium  gas  lasers 
aimed  at  photodetectors  across  the  impact  chamber  orthogonal  to 
and  intersecting  the  line  of  flight  of  the  projectile. 
Photomultipliers  monitor  the  laser  outputs  through  circular 
apertures  and  narrow  band  filters.  When  light  interruption 
occurs  at  each  station  because  of  projectile  passage,  a  sharp 
change  of  voltage  level  is  converted  into  a  trigger  signal  of 
sufficient  voltage  to  initiate  a  Field  Emission  Corporation 
30  nsec  dual  flash  x-ray  unit.  The  x-ray  flash  exposes  a  Pola¬ 
roid  film  plate  on  the  opposite  side  of  the  chamber  by  means  of 
a  fluorescing  intensifier  screen. 

Spacing  between  the  two  x-ray  field  position  referenced  is 
indicated  by  fiducial  wires  which  are  measured  by  an  optical 
comparator  to  within  0.1  mm.  Measurements  of  the  impactor  face 
position  relative  to  the  window  fiducials  allow  determination 
of  the  projectile  position  and  distance  traveled  over  the  time 
interval  measured  between  flash  exposures.  Figure  7  is  an 
example  of  the  shadowgraphs  taken  at  the  two  x-ray  stations 
showing  the  projectile  in  free  flight  before  impact. 

Serving  as  a  backup  technique,  a  second  method  is  employed 
to  measure  impactor  velocity  by  recording  the  time  interval 
between  the  first  x-ray  flash  and  the  event  of  shorting  of 
a  rear  surface  shock  velocity  pin.  The  impactor  and  target 
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positions  are  measured  from  the  x-ray  shadowgraphs .  Variations 
in  measurements  between  the  two  techniques  are  usually  less 
than  0.05%. 


IMPAf TOR  TARGFT 


. . .  >i  ii  *  —  t  i  vm 

STATION  #1  STATION  #2 


Figure  7  X-Ray  Shadowgraphs  of  Projectile 
Before  Impact 


The  early  experimentation  was  somewhat  different  from  that 
described.  A  single  laser  system  was  the  x-ray  trigger  source, 
and  the  x-ray  flashes  then  occured  at  preset  delay  times. 

During  low  and  intermediate  launch  velocity  tests,  the  system 
was  'v  95%  reliable.  But  in  the  highest  velocity  (7-8  km/sec) 
tests,  the  x-rays  were  nearly  always  pre triggered,  resulting  in 
the  x-ray  flashes  occurring  before  the  arrival  of  the  projectile. 
The  hugoniot  data  recorded  -  the  shock  wave  transit  times  - 
are  not  totally  invalidated,  when  the  impact  velocity  is  not 
recorded  directly.  Through  ballistics  the  light  gas  gun 
muzzle  velocity  is  well  documented  and  cross  correlated 
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with  secondary  parameters.  It  is  posible  through  the  use 
of  the  leading  data  and  related  cross  curves  of  hydrogen 
pressure,  powder  load,  and  piston  velocity  history,  to 
estimate  the  impact  velocity  to  within  +2%. 

The  target  is  located  approximately  60  cm  from  the  launch  tube 
muzzle  and  is  included  in  the  #2  x-ray  field  of  view.  Record¬ 
ings  of  the  shock  wave  transit  time  in  the  target  are  made 

using  four  coaxial  self-shorting  pins  as  sensors.  Shorting 
of  a  pin  results  in  a  sharply  rising  current  to  ground  which 
produces  a  signal  across  the  time-interval-meter  input  termina¬ 
tion  resistors.  The  pin  supply  circuit  is  so  designed  that 
each  pin  signal  can  be  seen  on  three  output  lines ,  free  of  any 
reflections  or  ringing  for  several  hundred  nanoseconds.  The 
Individual  circuits  are  "tuned"  by  the  use  of  trimmer  capaci¬ 
tors  so  that  the  risetime  of  each  signal  is  1.0  ±  0.1  nsec  to 
12  volts.  The  time-interval-meter  trigger  level  is  set  pre¬ 
cisely  at  1  volt.  Thus  the  combined  mechanical-electronic 
signal  system  can  make  use  of  the  ±  1/2  nsec  resolution  of  the 
time  recording  instruments. 

The  shock  wave  transit  time-interval-meters  are  Eldorado 
Models  793.  These  counters,  which  have  a  specified  time 
resolution  of  ±  1/2  nsec,  may  be  read  digitally  to  the  nearest 
nanosecond.  They  require  an  input  signal  of  1  volt  minimum  and 
a  risetime  of  approximately  1  nsec.  Although  instrument 

4 

stability  is  specified  to  be  one  part  in  10  for  long  term  use 

6 

24  hours)  and  five  parts  in  10  for  short  term  ('v  minutes)  , 
in  actual  practice  the  instruments  are  calibrated  prior  to 
each  shot  over  a  period  of  about  30  minutes.  The  shot  is  then 
fired  within  five  minutes  of  completion  of  the  calibration 
procedure . 
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The  coaxial  self-shorting  pins  employed  as  sensors  consist  of 
an  outer  sheath  of  brass  surrounding  a  teflon  sleeve  and  a  copper 
inner  conductor.  The  pins ,  connected  to  RG-174  50  Q  cable 
by  soldered  joints,  are  made  self-shorting  by  the  placement  of 
a  brass  cap  over  the  sensing  end  which  leaves  a  small  gap  (on 
the  order  of  0.050  ±  0.015  mm)  between  its  inner  face  and  the 
flat  end  of  the  inner  conductor.  Pin  gaps  are  recorded  by 
x-ray  shadowgraphy  of  which  Figure  8  is  an  example,  and  the 
measurements  taken  from  the  film  are  employed  in  the  shock 
velocity  calculations  for  corrections  as  pin  closure  times . 


Figure  8  X-Ray  Shadowgraph  of  Coaxial  Self-Shorting  Pin 


Model  Ca-1039,  Edgerton,  Germeshausen  and  Grier,  Santa 
Barbara,  California 
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Careful  attention  is  given  to  the  electrical  characteristics 
of  the  coaxial  self  shorting  pins.  A  Tektronix  1S2  time-domain- 
ref  lectome  ter  is  used  to  detect  slight  changes  in  the  reactance 
of  the  pins,  and  serve  as  an  indicator  of  probable  pin 
failure  in  the  vacuum  environment. 

Impactor  striking  tilt  is,  in  the  present  system,  recorded 
with  the  Eldorado  time-interval-meters.  Most  of  the  experiments 
reported  here  employed  a  Tektronix  519  oscilloscope  as  the 
impact  tilt  recorder.  The  20  nsec/' cm  standard  sweep  rate  was 
found  to  be  nonlinear  over  the  first  centimeter  of  the  CRT 
graticule.  Calibration  plots  were  periodically  generated 
and  employed  to  correct  the  recordings . 

Planarity  of  the  shock  wave  induced  in  the  target  is  depend¬ 
ent  on  the  impactor  flatness  at  impact.  The  impactor  surface  is 

.  -3 

machine  lapped  and  then  hand  polished  flat  to  0.5  x  10  mm. 

Tests  of  the  effect  of  launching  on  surface  flatness  were 

performed  by  impacting  into  mirrors  to  optically  observe  the 

impact  linearity.  With  impactors  of  6061-T6  aluminum,  the 

surface  curvature  over  a  44.5  mm  diameter  (length/diameter  - 

1:45)  at  impact  is  less  than  the  5  nsec  time  resolution  of  the 

optical  recorder  at  a  launch  velocity  of  ^  2  km/sec.  Similar 

tests  with  21  mm  diameter  impactors  at  launch  velocities  of  'v 

6  mm/sec  also  showed  <  5  nsec  curvature.  The  impactor  diameter 

for  equation  of  state  experiments  is  nominally  21  mm  and  the 

L/D  'v  1:12.  The  expected  launch  induced  curvature  for  the 

copper  and  Fansteel-77  impactors  is  <  1  nsec  because  of  the 

greater  stiffness  and  the  larger  L/D. 

Impactor  tilt  relative  to  the  target  specimen  front  surface 
is  sensitive  to  launch  tube  linearity  and  sabot  alignment  as 
well  as  to  target  alignment,  even  though  the  projectile  flies 
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free  from  the  launch  tube  before  impact.  The  capability  to 
adjust  the  target  position  and  perpendicularity  relative  to 
the  launch  tube  centerline  by  an  optical  technique  brings  the 
average  tilt  at  impact  to  approximately  0.05  radians  (approxi¬ 
mately  20  nsec  of  tilt  at  an  impact  velocity  of  6  km/sec) . 

Because  of  the  comparatively  gentle  acceleration  of  the  pro¬ 
jectile  to  its  terminal  velocity,  the  impactor  plate  is  not 
shock  stressed  or  shock  heated.  In  addition,  free  flight  in 
an  evacuated  range  precludes  aerodynamic  heating.  The  estimated 
temperature  rise  during  launch  is  of  the  order  of  1°C.  This 
temperature  increase  does  not  significantly  affect  the  flatness 
of  the  impactor. 


Target  Design  and  Construction 

The  photograph  and  drawing  of  Figure  9  illustrates  the  basic  fea¬ 
tures  of  the  target  design.  Two  coaxial  self-shorting  pir.r  (sta¬ 
tions  1  and  A)  are  passed  by  the  edge  of  the  specimen  with  their 
cap  faces  exactly  in  the  plane  of  the  specimen  impact  surface. 
These  pins  are  used  to  initiate  the  timing  for  the  shock  wave 
transit  and  to  measure  impactor  tilt  as  the  time  interval 
differences  between  their  respective  closures.  Two  rear 
surface  pins ,  stations  B  and  C,  (or  B  only  depending  upon  the 
impactor  diameter)  are  mounted  in  line  with  the  tilt  pins  to 
record  the  shock  wave  arrival  at  the  rear  surface.  All  four 
pins  are  mounted  in  a  guide  fixture  which  assures  accurate 
geometrical  spacing.  The  tilt  pins  are  fixed  in  position  by 
dimensionally  stable  epoxy,  while  the  rear  surface  pins  are 
spring  loaded  against  the  target.  The  pin  retainer  and  cable 
bracket  lend  rigidity  to  the  assembly  preventing  accidental 
damage  to  the  pin  shafts. 
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Two  target  sizes  were  chosen,  each  to  be  used  over  a  given 
stress  range.  The  four-pin  target  required  a  21  mm  diameter 
impactor  which  was  limited  to  a  maximum  impact  velocity  of 
about  6.5  km/sec.  Single-rear-surface-pin  targets  requiring 

II. 7  mm  diameter  impactors  were  used  for  launch  velocities  up 
to  %  8  km/sec.  The  four-pin  targets,  in  conjunction  with 
four  Eldorado  one-nsec  time  interval  meters,  allow  the 
determination  of  two  sets  of  measurements,  yielding  four 
values  of  shock  wave  transit  time.  As  described  in  Section 

III,  these  four  values  are  combined  to  calculate  a  single 
shock  wave  velocity  in  the  specimen,  and  by  a  system  of  cross¬ 
checking,  to  provide  an  indication  of  the  accuracy  of  the 
measurement. 

Throughout  the  program,  the  edge  rarefaction  angle  (a  in  Figure 
9b)  was  assumed  to  remain  constant  at  tan  a,  -  1.0  which  is  in 
agreement  with  a  Poisson' a  ratio  of  ^  0.3,  typical  for  rocks 
under  hydrostatic  stress.  Similarly,  the  target-impactor 
thickness  ratios  of  3.3  and  3.8  for  the  two  target  sizes  were 
established  so  as  to  be  thinner  than  the  point  at  which  the 
rear  surface  (forward  going)  elastic  rarefaction  overtakes  the 
shock  front. 

Except  for  the  granite  pnd  most  of  tuff  targets,  the  specimens 
were  machined  into  solid  rectangles  where  the  dimension  along 
the  pin  array  was  17.6  mm  and  9.7  mm  for  the  two  target  sizes. 
The  base  line  (X3  in  Figure  9b)  was  19.2  mm  and  11.2  mm  for  the 
intermediate  and  the  high  velocity  targets,  respectively. 

The  impact  and  rear  surfaces  were  machine  lapped  and,  where 
possible,  hand  polished  to  as  fine  a  finish  as  the  material 
would  take.  The  lapping  procedure  also  served  to  enhance  the 
surface  flatness  and  parallelism.  These  specimen  surfaces  were 
generally  flat  and  parallel  to  within  5*10-3  mm. 
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Target  specimen  thicknesses  at  the  pin  stations  were  mea¬ 
sured  with  a  Zeiss  light-section  microscope  employed  as  a  com¬ 
parator.  The  references  were  precision  laboratory-grade  gage 
blocks.  Use  of  the  light-section  microscope  avoids  the 
problem  of  an  indicator  marring  the  specimen  surfaces . 
Specimens  were  typically  v  5.8  mm  and  3.9  mm  thick. 


Tuff  Targets 

The  tuff  specimens  presented  several  special  problems ,  since  the 
samples  are  ^  20%  porous  and  the  rock  is  hygroscopic.  The  mea¬ 
surements,  construction,  and  storage  were  performed  at  atmospheric 
conditions  after  having  determined  that  short  5  min)  exposure 
to  vacuum  would  almost  completely  dry  the  rock.  Exposure  to 
vacuum  for  at  least  a  10  min  duration  always  proceeded  the  shots. 

The  requirement  later  in  the  program  to  simulate  in  situ  conditions 
made  it  necessary  to  maintain  a  liquid  environment  about  the  porous 
tuff  specimen  for  a  period  long  enough  to  ensure  saturation.  To 
insure  complete  saturation,  the  specimen  had  to  be  submerged  in 
water  for  approximately  48  hours  (i.e.  at  least  ^  95%  of  the 
water  was  absorbed  within  48  hours;  the  final  'u  5%  took  up  to 
2  weeks  at  atmospheric  pressure) .  For  the  experiments  the 
liquid  enclosure  had  to  also  withstand  an  external  vacuum  in 
the  target  chamber.  The  basic  target  design  will  be  described 
later  as  the  stepped  model  with  the  addition  of  a  cell  cover 
plate  and  appropriate  water  seals. 

After  the  target  was  epoxied  in  the  mounting  block,  a  ^  0.1  mm 
layer  of  epoxy  was  cast  over  the  entire  front  surface  of  the  tar¬ 
get  specimen  and  mounting  block.  The  casting  was  made  using  shim- 
stock  spacers  between  the  casting  fixture  and  the  target  block 
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surfaces.  Thickness  measurements  made  with  the  Zeiss  liqht  sec- 

tion  microscope  showed  that  the  resulting  epoxy  surface  was  with- 

-3 

in  10  mm  of  being  parallel  with  the  target  surface.  Microscopic 
examination  of  epoxy  coated  specimens  showed  epoxy  penetration 
to  be  less  than  0.03  mm.  The  effect  upon  the  velocity  of  a  shock 
wave,  initiated  through  the  shim  was  shown  by  a  computer  hydrody¬ 
namic  calculation  to  be  less  than  a  0.1%  retardation. 

Target  specimen  saturation  was  then  begun  by  evacuating  the 
assembly  and  introducing  deaireated  water  into  the  evacuated 
specimen  cell.  Without  this  procedure,  capillary  action  forced 
trapped  air  toward  the  thin  epoxy  front  surface  shim  causing 
the  shim  to  release  from  the  tuff  specimen  impact  surface. 

The  rear  surface  cell  cover  was  epoxied  in  place  with  the  tar¬ 
get  specimen  in  the  wet  condition.  Fill  tubes  were  provided 
on  the  cell  through  the  cover  plate. 


Grout  Targets 

The  tuff  impedance-matching  grout,  V44D,  was  rather  plastic  at 
room  temperature  and  would  deform  about  the  pin  caps  under  the 
light  but  continuous  stress  of  the  pin  retaining  spring.  To 
limit  the  deformation  as  much  as  possible  and  maintain  the  tar¬ 
get  thickness  as  recorded  optically,  a  brass  shim  was  placed 
under  the  rear  surface  pin  caps  thus  distributing  the  pressure 
over  a  larger  area.  The  transit  time  measurements  were  corrected 
for  the  shim  as  described  in  Section  III. 


Alluvium  Targets 

The  alluvium  specimens  were  prepared  in  much  the  same  manner  as 
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described  by  Anderson,  G.  D.  and  A.  L.  Fahrenbruch ^ ,  except 
that  the  pressed  sample  surfaces  were  not  covered  with  aluminized 
mylar. 

To  avoid  the  effects  of  an  unknown  amount  of  water  being  absorbed 
from  the  air,  sample  preparation  was  carried  out  under  controlled 
atmospheric  conditions.  The  loose,  finely  particulate  material  was 
screened  (No.  40  mesh)  and  placed  overnight  in  a  drying  oven. 

Upon  removal  from  the  oven  all  samples  were  stored  in  a  glove 
box  set  up  to  maintain  an  atmosphere  of  less  than  1%  absolute 
humidity.  Specimens  were  exposed  to  room  air  for  no  longer 
than  5  minutes  thereafter. 

The  alluvium  was  pressed  into  a  steel  ring  with  inside  diameter 
of  15.2  mm  and  thickness  was  5.1  mm.  The  front  and  rear  surfaces 
of  the  ring  were  machined  and  lapped  flat  and  parallel  to  the  usu¬ 
al  tolerances,  and  each  specimen  ring  was  precisely  measured  and 
weighed.  The  alluvium  sample  was  then  compacted  into  the  ring 
at  the  desired  specimen  density.  Tests  conducted  on  several 
trial  specimens  showed  that  the  surface  of  the  compacted  material 
had  sufficient  rigidity  to  not  require  a  shim  under  the  rear  sur¬ 
face  pin  caps . 


The  front  surface  pins  were  placed  in  holes  drilled  through  the 
ring  so  that  the  pins  are  aligned  in  a  manner  similar  to  that 
in  the  regular  target  construction.  Finally,  the  ring  assembly 
was  epoxied  to  a  modified  target  mounting  block. 


Stepped  Target  Design 


Tests  conducted  during  the  early  part  of  the  research  program 
used  a  different  target  design.  Figure  10  illustrates  this  de- 
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sign.  Originally,  the  tilt  pins  (stations  1  and  A)  were  not  in 
the  plane  of  the  impact  surface;  rather,  they  rested  on  two  wings 
of  the  target  of  ^  0.76  mm  thickness.  The  B  and  C  pins  were 
placed  on  steps  of  thickness  3.30  and  5.8  mm  respectively.  The 
pin  guide  fixture  was  about  4-5  mm  above  the  target  C  step. 


BRASS  MOUNTING  BLOCK 


COAXIAL  SHORTING 
PINS 


PLEXIGLAS 
PIN  GUIDE 


Figure  10  Early  Equation  of  State  Target  Design 
(Step  Target) 

Shock  wave  velocities  were  measured  through  two  different  thick¬ 
nesses  of  material.  Thus,  non-steady  state  wave  behavior  in 
the  material  could  be  measured,  if  present.  This  more  compli¬ 
cated  and  time  consuming  target  system  was  discontinued  in  favor 
of  the  system  previously  described  when,  within  the  resolution 
of  the  experiment,  no  significant  differences  in  wave  velocity 
were  measured. 
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Machining  of  the  stepped  target  configuration  was  difficult, 
and  pin  placement  with  the  raised  pin  guide  could  allow  the 
pin  to  slant  or  be  displaced  from  the  perpendicular  to  the  tar¬ 
get  surface.  In  addition,  higher  tilt  measuring  accuracy  could 
be  gained  by  removing  possible  specimen  material  effects  in  the 
wings  and  increasing  by  about  5  mm  the  tilt  measurement  base 
line . 
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SECTION  II 


DATA  ANALYSIS 


Analysis  of  the  experimental  data  obtained  for  the  program 
was  based  upon  the  impedance  match  solution  for  the  determina¬ 
tion  of  particle  velocity  and  pressure  of  the  shock  state 
achieved  in  the  specimen  during  each  test.  The  requisite 
information  in  this  analysis  is  the  shock  wave  velocity,  the 
impact  velocity,  and  the  linear  form  of  the  hugoniot  of  the 
impactor  material.  A  detailed  analytical  discussion  of  the 
approach  is  presented  in  the  Appendix . 


Shock  Wave  Velocity 

The  measurements  relevant  to  shock  wave  velocity  are  the  tar¬ 
get  thickness  and  the  shock  wave  transit  time  interval.  To 
determine  as  precisely  as  possible  the  shock  wave  transit  time, 
refinements  may  be  applied  to  the  recorded  time  interval. 

Two  sources  of  refinement  to  the  measurement  are: 

(i)  Inclusion  of  the  effects  of  shock  wave  tilt 

resulting  from  non-coplanar  impact  of  the  target 
and  impactor  surfaces. 

(ii)  Correction  for  the  differences  of  closing  times 

of  coaxial  pins  with  varying  gaps,  due  to  impactor- 
pin  and  target-pin  interactions. 
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The  coaxial  shorting  pin  design  includes  a  cap  over  the  impact 
end  of  the  pin  which  has  a  slightly  recessed  rear  surface.  The 
gap  between  the  cap  and  the  inner  electrode  is  electrically 
insulating.  When  a  large  amplitude  stress  wave  reaches  the 
cap  face,  the  cap  is  set  into  motion  and  the  gap  is  assumed 
to  be  closed  at  a  rate  equal  to  the  free-surface  velocity  of 
the  cap  material.  Correction  for  shock  transit  through  a  rear 
surface  shim  is  made  by  simply  defining  the  pin  cap  thickness 
as  the  sum  of  the  shim  and  measured  pin  cap  thicknesses. 

The  system  of  analysis  has  been  incorporated  in  a  computer 
program  SHOVEL  -  presented  in  the  Appendix,  and  is  perma¬ 
nently  stored  in  a  General  Electric  Company  computer  time¬ 
sharing  service  library. 

Major  sources  of  experimental  error  include  unrecorded  impact 
velocity,  pin  gap  closure  velocities  and  the  impactor  tilt. 

An  unrecorded  impact  velocity  results  in  an  uncertainty  of  'v  4% 
in  the  specimen  particle  velocity  but  affects  the  shock  velocity 
measurement  insignificantly,  i.e.  <  0.1%.  Pin  cap  closure 
velocity  is  subject  to  the  errors  of  the  assumption  of 
upin  free  surface  =  2upin‘  In  Practi°e,  however,  this  assumption 

has  a  maximum  effect  on  Ug  of  ^  0.1%  since  the  pin  closure  time 
is  on  the  order  of  100  times  smaller  than  the  shock  transit  time. 


The  effect  of  impactor  tilt  on  U  can  be  quite  large  as  seen  by 

inspection  of  Equation  A-18.  However  a  measurement  error  of  100% 

results  an  uncertainty  of  U  of  only  ^  2%.  The  system  employed 

s 

for  tilt  measurement  for  these  experiments  is  accurate  to  within 
0.0015  radians.  The  combined  effect  of  these  errors  results  in 
an  uncertainty  in  U  of  between  0.2  to  0.5%,  depending  upon  the 

o 

target  diameter  and  the  shock  wave  transit  time. 
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As  has  been  shown,  two  values  of  shock  velocity  and  hence  two  sets 
of  values  of  particle  velocity,  pressure  and  specific  volume  are 
determined  with  each  target  in  the  four-pin  configuration.  The 
data  presented  in  Section  IV  are  shown  as  averages  of  the  shock 
velocity  and  particle  velocity  between  the  B  and  C  station  mea¬ 
surements  for  each  test  whenever  possible.  The  "internal"  (or 
individual  test)  data  scatter  is  listed  and  displayed  graphically 
as  error  bars  on  all  shots  where  averaging  was  possible.  The  in¬ 
ternal  scatter  in  most  instances  is  outside  the  experimental  sys¬ 
tematic  errors  but  in  those  cases  where  partial  system  failure 
occurred  the  error  bars  include  the  system  error. 

The  linear  form  U  (u  )  of  the  least  squares  analysis  was 

s  P 

employed  to  represent  analytically  the  hugoniots  of  the  materials 
tested  in  this  program.  Where  the  data  suggest  systematic  devia¬ 
tions  from  linearity,  the  analytical  form  was  considered  to  be 
represented  by  additional  linear  segments .  Thus  the  reported 
hugoniot  equations  are  followed  by  a  designation  of  the  range 
over  which  the  least  squares  analysis  was  carried  out.  The 
standard  errors  (S.E.)  are  also  listed  to  indicated  the  least 
mean  square  of  the  deviations  of  the  dependent  variable  (U  ) 
from  the  computed  linear  fit. 
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SECTION  III 


EXPERIMENTAL  RESULTS 


Presented  here  are  the  results  of  the  experiments  and  data  an¬ 
alysis  performed  during  the  program.  Although  detailed  inter¬ 
pretation  of  the  data  is  beyond  the  scope  of  this  work,  where 
possible,  comparisons  of  the  results  of  other  researchers  which 
are  considered  to  be  of  interest  are  made. 


Impactor  Materials  Hugoniot  Data 

The  hugoniot  data  for  Fansteel-77  are  listed  in  Table  I. 

Indicated  in  the  table  are  those  shots  performed  under  General 
Motors,  NASA,  and  BRL  sponsorship.  The  shot  numbering  indicates 
the  laboratory  at  which  the  tests  were  performed:  C,  at  GM 
Defense  Research  Laboratories,  S,  at  Manufacturing  Development, 

GM  Technical  Center.  The  seventeen  shots  listed  establish  with 
precision  the  Fansteel-77  hugoniot  over  the  pressure  range  attain 
able  in  the  rock  experiments . 

The  physical  properties  of  Fansteel-77,  which  is  a  powder  metal¬ 
lurgy  product,  are  presented  in  Table  II.  Microscopic  examina¬ 
tion  reveals  that  the  tungsten  component  is  granular  with  the 
average  grain  diameter  ^  1  y .  X-ray  diffraction  measurements 
show  that  the  intragranular  matrix  is  an  alloy  of  'v  70%  nickel 
with  ^  30%  copper. 
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TABLE  II 


CHEMICAL  AND  PHYSICAL  PROPERTIES 
OF  FANSTEEL-77 


Chemical  Properties 


Element 

W 

Cu 

Ni 


Wt  % 

90  ±  1 
3.8  ±  0.6 
6.1  ±  0.2 


Physical  Properties 


Ultimate  Compressive  Strength 
Ultimate  Tensile  Strength 
%  Elongation 

Yield  Strength  (0.2%  elong,') 
Density 

Sonic  Velocity*  Longitudinal 
Sonic  Velocity*  Shear 


176,000  psi 
98,000  psi(min) 
2%  (min) 
85,000  psi(min) 
17.01  gm/cm^ 
5.049  km/sec 
2.765  km/sec 


Calculated  Parameters 


Bulk  Sound  Velocity 
Bulk  Modulus 
Bar  Wave  Velocity 
Young's  Modulus 
Poisson's  Ratio 
Rigidity  Modulus 
Lame's  Parameter 
Rayleigh  Wave  Velocity 


3.911  km/sec 
2.602  Mb 
4.432  km/sec 
3.341  Mb 
0.286 
2.340  Mb 
1.042  Mb 
2.628  km/sec 
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The  sonic  velocity  measurements  were  performed  at  Manufacturing 
Development  by  the  pulse  superposition  technique .  Using 
the  manufacturers  value  of  Poisson's  ratio  the  bulk  sound 
speed  is  3.88  km/sec.  The  Poisson's  ratio  needed  to  bring  into 
agreement  the  linear  Ug-u^  hugoniot  intercept  of  4.00  ^  km/sec 
for  the  present  alloy  is  0.31. 

The  hugoniot  data  for  Fansteel-77  (with  a  measured  density  of 

3 

17.01  ±  0.01  gm/cm  )  in  the  U  -u  plane  may  be  described  by 

s  p 

the  linear  equation: 

U  =  4.008  +  1.262  u  km/sec,  0.4  <  u  <3.5  km/sec 
s  P  P 

S.E.(U  )  =  0.036  km/sec 
s 

Tests  exhibiting  less  than  1%  internal  (B  vs  C  Station)  scatter 

in  U  were  employed  in  the  least  squares  analysis  (7  tests) . 
s 

The  effect  of  employing  the  additional  twelve  shots  changes 
the  slope  by  only  0.03%.  The  data  are  displayed  along  with 
the  OFHC  copper  data  in  Figures  11  and  12. 

OFHC  copper  (99.99%  purity)  was  employed  as  the  low  and  inter¬ 
mediate  pressure  standard  (i.e.  ^  0.1  to  ^  1.0  Mb)  in  the  hug¬ 
oniot  experiments  on  the  rocks  in  the  program.  A  number  of 
other  laboratories'  '  ’  have  performed  hugoniot  te  ts  for 
copper  which  are  in  generally  quite  good  agreement  with  this 
laboratory.  the  hugoniot  equation  for  the  present  data  in  the 
Ug-Up  plane  is  given  by: 

U  =  3.960  +  1.464  u  km/sec,  1.7  <  u  <4.7  km/sec 
s  P  P 

S.E.(Ug)  =  0.018  km/sec 
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All  shots  listed  exhibited  less  than  2%  B  to  C  station  scatter 

3 

in  U  .  The  standard  density  employed  was  8.930  gm/cm  .  Fig- 
s 

ure  12  illustrates  that  the  shock  hugoniot  in  pressure  of  copper 
is  slightly  more  than  half  that  of  Fansteel-77.  The  hugoniot 
data  are  presented  in  Table  III. 

The  third  standard  employed  to  achieve  the  lowest  pressure  data 
point  in  the  saturated  tuff  experiments  was  6061-T6  aluminum. 

The  hugoniot  of  Al'tshuler'  '  was  employed. 

U  =  5.380  +  1.350  u  km/sec 
s  p 

3 

The  impactor  density  was  2.72  gm/cm  . 

Climax  Stock  Granodiorite  (Granite) 

The  physical  and  petrographic  description  of  the  Climax  Stock 
Grandodiorite,  "Granite",  as  well  as  the  source  of  the  samples 

(9) 

is  presented  by  Barnes.  A  small  sample  of  the  original 

* 

granite  was  tested  independently,  and  showed  a  dry  effective 
porosity  of  1.75%.  The  measured  density  for  the  target  speci- 
mens  for  the  hugoniot  tests  was  2.65  ±  0.01  gm/cm  .  The  cal- 
culated  bulk  sound  speed  from  reference  '  is  in  the  range  of 
4.07  to  4.43  km/sec. 

Granite  was  the  first  rock  tested  in  the  program,  and  the 
target  design  used  was  the  stepped  configuration. 

The  data  presented  in  Table  IV  represent  those  tests  which  meet 
the  requirements  of  low  shock  tilt  and  internal  scatter.  Fig¬ 
ures  13  and  14  display  the  data  in  the  Ug-Up,  an<3  P-V  planes 

* 

Hornkohl  Laboratories ,  Inc.,  Bakersfield,  California 


33 


MANUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 

MSL-68-15 


x 

CO 


o 

CN 


>»co 

o 

vo 

(N 

CO 

a> 

o> 

CO 

CO 

rH 

o 

P  0 

in 

in 

vo 

O 

OV 

in 

CO 

n* 

r- 

rH 

o 

VO 

rH 

o 

r- 

rH 

in 

in 

in 

o 

O 

tr 

a  «  \ 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

£  £ 

CN 

co 

CO 

co 

n* 

in 

in 

in 

in 

<u  & 

Q  — 

r-H 

1 — 1 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

• — i 

rH 

o 

GO 

a\ 


II  II 

0  O 

>  Q- 


a 

a 


w 

Pi 

Pi 

o 

u 

u 

as 

En 

o 

« 

s 

g 

< 

Q 

Eh 

O 

H 

53 

O 

O 

D 

as 


1 _ 


Pi 


<U  >,— 
cH  p  a 

O  -H  <U 
Di-H  O  (0 

3  P  O  *> 

U  r-H  g 

to  at  .* 
Pi  >  —' 


P  r-H 
O  nJ 
P  -H 
O  P 
rt  0) 
Qi-P 
6  « 
H  £ 


1 1 

CO 

VO 

co 

o 

l'' 

in 

VO 

CO 

vo 

Ch 

in 

3  x 

CM 

VO 

n* 

CO 

o 

00 

00 

00 

vo 

vo 

co 

>  Hn 

o  g 
>  o 

CO 

r- 

r- 

r- 

VO 

vo 

VO 

VO 

vo 

VO 

vo 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

•H 

o 

o 

co 

C" 

a> 

N* 

rH 

CD 

o\ 

vo 

vo 

> 

p  3 

in 

co 

rH 

CO 

co 

CM 

CN 

CN 

o> 

vo 

N 

(TJ  rH 

co 

00 

vo 

uo 

CN 

rH 

rH 

rH 

<T> 

<J\ 

r- 

vo 

> 

H  O 

r- 

vo 

vo 

vo 

VO 

VO 

VO 

VO 

in 

in 

in 

in 

0)  > 

• 

■ 

• 

• 

• 

• 

• 

r 

• 

» 

Pi 

o 

o 

o 

o 

o 

O 

O 

O 

o 

o 

o 

o 

<T> 

•H 

VO 

CM 

O 

co 

o 

in 

vo 

CO 

CM 

CO 

UO 

CO 

CN 

OV 

o 

CN 

00 

in 

o 

CT\ 

in 

00 

O 

00 

<7\ 

cn 

TT 

o 

in 

O 

IH 

rH 

CM 

CM 

CN 

CM 

CN 

CO 

co  • 

TT 

N* 

CO 

co 

rH 

00 

r- 

0^ 

CTV 

rH 

00 

rH 

00 

N* 

rH 

■vr 

co 

O 

rH 

O 

rr 

VO 

00 

r** 

VO 

r- 

co 

VO 

00 

CN 

in 

in 

m 

av 

<T\ 

CO 

vo 

rH 

CM 

CN 

CN 

CO 

CO 

CO 

co 

co 

CO 

N* 

P  u 

X  -p  at 

00 

in 

CM 

CTi 

vo 

cr» 

VO 

n 

CN 

o 

m 

U  o  ui 

vo 

rH 

rH 

vo 

VO 

O'. 

00 

o 

ov 

CD 

in  O  O  \ 

N* 

CO 

o 

VO 

o 

Is- 

CO 

co 

r- 

1  #C  *H  g 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

V)  i)  X 

vo 

r** 

00 

00 

o\ 

CTi 

av 

o\ 

cr> 

o 

o 

>  — 

rH 

rH 

P  P  o 

OHO) 

VO 

uo 

<Tt 

co 

00 

rH 

CD 

in 

CN 

lit  U  to 

CM 

r- 

CN 

CO 

rH 

cn 

CM 

CM 

CO 

o 

P-V 

J>  ft  o  \ 

"<r 

00 

crv 

vo 

o 

O 

rH 

<T\ 

CT\ 

g  i-h  e 

■ 

• 

• 

• 

• 

• 

« 

• 

• 

• 

* 

h  at  32 
> 

CO 

co 

r-* 

vo 

r- 

vo 

r- 

3WiH(/33P33W3WW 

Uti<PiOOaOfciU(ii(ii 


jj 

rH 

O 

CO 

in 

CM 

vo 

CO 

o 

a\ 

CM 

OV 

rH 

O 

CO 

o\ 

rH 

00 

vo 

a> 

N* 

00 

rH 

VO 

o 

CN 

1 

1 

o 

i 

av 

o 

CM 

1 

| 

fH 

CO 

rH 

1 

O 

rH 

1 

o 

w 

C/3 

rH 

1 

o 

w 

i 

o 

• — 1 

I 

O 

rH 

I 

O 

CO 

cn 

1 

cn 

34 


Tests  performed  under  contract  NAS2-3427 


TABLE  IV 

HUGONIOT  DATA  FOR  CLIMAX  STOCK  GRANODIORITE  (GRANITE) 
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Single  pin  data 


SHOCK  VELOCITY  -  KM/SEC 
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PARTICLE  VELOCITY  -  KM/SEC 


Figure  13 


Shock  Velocity  vs  Particle  Velocity  for 
Climax  Stock  Granodiorite 


PRESSURE 

O-0CL  !•( 


MATERIAL  - 


GRANDDIORITE 


SPECIFIC  VOLUME  -  CC/GM 

Figure  14  Pressure  vs  Specific  Volume  for  Climax 
Stock  Granodiorite 
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respectively.  Points  labeled  single  ;.in  values  are  those 
shots  employing  only  a  single  shock  trai.it  time  pin. 

Error  bars  in  the  u^  axis  represent  the  extremes  of  possible 
particle  velocities  induced  by  the  impactor  whose  velocity  was 
determined  by  measured  gun  firing  parameters. 

The  linear  hugoniot  equation  is  given  by: 

U  =  2.169  +  1.644  u  km/sec,  3.1  <  u  <6.0  km/sec 
s  P  P 

S,E.(Us)  =  0.190  km/sec  for  10  points. 

Figure  15  is  a  comparison  of  the  present  results  with  the  data 
of  McQueen^^  (for  Westerly,  R.I.  granite)  and  the  LRL  Compiler 
for  several  NTS  Area  15  granites.  The  hugoniot  equation  for  the 
high  density  portion  of  McQueen's  Westerly  granite  is: 

Us  =  2.102  +  1.629  Up  km/sec,  2.5  <  u^  <  4.1  km/sec 

S.E.(U  )  =  0.052  km/sec  for  8  points 

fa 

The  least  squares  linear  equation  for  the  LRL  Compiler  data  is 
given  by: 

Ug  =  2.594  +  1.476  u  km/sec,  2.5  <  u  <3.9  km/sec 

P  P 

S.E.(Ug)  =  0.056  km/sec  for  4  data  points. 

The  large  standard  error  of  the  GM  data  for  this  material  is 
attributable  to  the  many  experimental  problems  encountered 
~arly  in  the  program.  The  fact  that  the  hugoniots  of  both  LRL 
and  McQueen  compare  vei  closely  with  the  GM  results,  even  for 
different  granites,  suggests  the  larger  error  does  not  signifi¬ 
cantly  effect  the  results. 
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A  collection  of  all  the  granite  high  pressure  phase  data  yields 
the  following  equation: 

U  =  1.994  +  1.673  u  km/sec,  2.5  <  u  <  6.0  km/sec 
s  p  p 

S.E.(U  )  =  0.149  km/sec  for  22  points. 

This  combining  of  data  appears  reasonable  for  a  general  formula¬ 
tion  of  the  high  pressure  "granite"  hugoniot  based  on  the  initial 
densities  and  the  mineral  compositions  shown  below. 


GM 

McQueen  ^ 

LRL  Compiler 

Source 

NTS  Area  15 

Westerly,  R.I. 

NTS  Area  15 

.  3 

Ave.  Density (gm/ cm  ) 

2.65 

2.628 

2.677 

Composition  (Wt  %) 

Quartz 

21-22% 

27.5% 

19.5% 

Plagioclase 

46-49% 

31.4% 

26.5% 

Orthoclase 

17-20% 

35.4% 

47.0% 

Total  Feldspar 

63-69% 

66.3% 

73.5% 

Although  there  ware  no  tests  performed  in  the  lower  phase  regime, 
the  following  are  the  hugoniots  obtained  by  McQueen  and  LRL  re¬ 
spectively  for  this  region: 

U  =  4.926  +  0.371  u  km/sec,  1.0  <  u  <2.2  km/sec 
“  P  P 

S.E.(U  )  =  0.019  km/sec  for  11  points. 

O 

U  =  5.354  +  0.159  u^  km/sec,  0.4  <  u  <2.2  km/sec 
is  P  p 

S.E.(U  )  =  0.134  km/sec  for  9  points. 


(11) 
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Only  the  higher  density  of  the  two  sets  of  LRL  Compiler 
granite  data  are  included  here  and  in  the  figures.  The 
granites  compared  here  appear  to  exhibit  similar  hugoniots 
above  the  high  density  phase  transition  but  appear  to  differ 
somewhat  more  below  the  transition. 


Rainier  Mesa  Tuff  (Vacuum  dry) 

The  physical  and  petrographic  descriptions  for  Rainier  Mesa 

(9 ) 

tuff  are  reported  by  Barnes.  The  source  of  the  samples 

is  also  given  in  this  description. 

The  target  configuration  for  the  dry  tuff  and  most  of  the  water 
saturated  tuff  experiments  was  the  stepped  design.  It  was 
necessary  to  select  target  specimens  having  good  uniformity 
and  small  granularity  in  the  tilt  pin  steps. 

Hugoniot  data  for  the  tuff  initially  in  the  vacuum  dry 
condition  are  presented  in  Table  V,  and  Figures  16  and  17. 

The  linear  form  of  the  hugoniot  is  represented  by  the  least 
squares  equation: 

U  =  1.049  +  1.528  u  km/sec,  2.3  <  u  <  6.5  km/sec 
s  p  p 

S.E.(U  )  =  0.086  km/sec  for  8  data  points, 
s 

3 

The  initial  density  was  1.76  ±  0.02  gm/cm  as  measured  prior  to 
the  tests. 

Comparison  of  the  present  data  with  selected  results  of  other 
researchers  is  illustrated  in  Figure  18.  Even  though  the 
specimens  of  Lombard,  Wiedermann ,  ^ ^ ^  Anderson,  and 


41 


\ 

i 


Single  Pin  Value 


PRESSURE 


MANUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 


MSL-68-15 


Bass  are  reported  to  be  from  NTS  Area  12,  the  Rainier  Mesa 

area,  there  is  considerable  variability  in  the  tuff's  physical 

properties.  The  density  and  porosity  of  the  comparison  data 

3 

were  chosen  to  be  in  the  ranges  pQ  ^  1.61  to  1.80  gm/cm  and 
^  20-30%  respectively. 

Nevertheless  the  data  scatter  is  small  as  the  figure  displays 
and  establishes  a  fairly  well  defined  hugoniot  when  taken 

together . 


A  least  squares  fit  to  the  selected  and  combined  data  of  the 
low  pressure  regime  is  defined  by  the  linear  equation: 


U  -  1.371  +  1.431  u  km/sec,  0.7  <  u  <2.5  km/sec 
s  P  P 

S.E.(U  )  =  0.116  km/sec  for  32  points, 
s 


There  is  a  slight  disconinuity  between  the  comparison  data  and 
the  present  linear  fit  to  the  higher  pressure  data.  This  break 
in  the  hugoniot  shows  up  best  in  the  P-V  plane.  Assuming  that 
the  data  shown  in  Figure  18  (see  also  Figure  19)  for  the  lower 
pressures  are  reliable  and  representative,  a  break  appears  to 

3 

occur  at  ^  0.190  Mb  and  specific  volume  of  0.30  cm  /gm.  Using 

3 

the  average  value  of  density  (PoAve  =  1-647  gm/cm  )  and  the 
slope  of  the  linear  fit  to  all  the  data,  the  following  equa¬ 
tion  appears  to  describe  a  zone  in  P  =  P(V)  which  represents 
the  range  of  the  hugoniot  states  of  the  various  tuffs  up  to  a 
pressure  of  ^  0.19  Mb. 


P 


CQ2  <VQ-V> _ 

V  -  1.431  (V  -V)l  2 

.  °  °  J 


46 


PRESSURE 


MANUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 


MSL-6S-15 

3 

where  V  =  0.607  cm  /gm,  and  Cq  may  range  from  0.132  to 
0.147  cm/ysec.  The  values  of  the  initial  densities  of  the 
various  tuffs  are  included  in  the  legend  of  Figure  19  . 

The  evidence  presented  hero  for  the  existence  of  the  phase 
change  in  the  dry  tuff  is  still  subject  to  dispute,  primarily 
because  of  the  question  of  the  reliability  of  the  data  in  such 
varied  polymineralic  rock  and  of  the  mechanism  which  causes 
the  spread  in  the  shock  velocity  measurements.  Unfortunately 
the  lowest  present  experimental  datum  is  either  just  above  or 
at  the  discontinuity  allowing  no  direct  comparisons  with  the 
lower  end  of  the  hugoniot.  The  comparison  made,  however,  be¬ 
tween  the  several  points  reported  by  Lomard,  Wiedermann,  and 
Bass  on  either  side  of  the  suggested  transition  point  appears 
to  be  in  agreement. 


Rainier  Mesa  Tuff  (Water  Saturated) 


The  design,  construction  and  testing  of  the  water  saturated 
tuff  targets  is  presented  in  Section  III  in  detail. 

Each  target  was  checked  for  water  leakage  in  vacuo  prior  to 
the  actual  hugoniot  experiment.  Several  pin  failures  due  to 
preshorting  did  occur;  these  are  indicated  in  Table  VI. 

Figures  20  and  21  graphically  present  the  data  for  saturated 

3 

tuff.  The  measured  initial  density  was  1.97  ±  0.02  gm/cm 
for  these  specimens.  The  data  are  best  represented  by  two 
line  segments.  The  data  are  amenable  to  description  by  a 
linear  segment  given  by: 
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U  =  1.755  +  1.575  u  km/sec,  3.0  i  u  <6.1  km/sec 
s  P  P 

S.E, (U  )  =  0.089  km/sec  for  6  points, 
s 

The  lower  segment  is  much  less  well  defined  from  the  present- 
data  alone,  but  the  linear  form  is  presented  for  convenience 
as : 


U  =  2.874  +  1.141  u  km/sec,  1.0  <  u  <  2.0  km/sec 
s  P  P 

S.E.(U  )  =  0.162  km/sec  for  5  points, 
s 

The  reasons  for  the  data  scatter  in  the  lower  region  are  pro¬ 
bably  several.  The  two  lowest  pressure  tests  experienced  a 
failure  due  to  water  preshorting  of  one  of  the  two  rear  sur¬ 
face  pins.  Also  at  low  shock  stress  and  shock  velocity,  i.e. 
velocities  below  the  sonic  velocities  of  the  basic  minerals, 
it  is  expected  that  inhomogeneities  of  the  specimens  will 
effect  wave  propagation,  as  will  specimen  to  specimen  variations 

Shock  wave  arrival  at  a  free  surface  in  porous  and  granular 
rock  materials  has  been  observed  to  be  irregular  and  non- 
planar.  With  such  few  data  available  the  representation 

cannot  be  considered  to  be  an  exact  description.  Based  upon 
the  obvious  scatter  of  the  data  of  other  authors  the  velocities 
recorded  are  probably  representative  of  the  bulk  material. 

Figures  22  and  23  are  comparison  plots  of  the  present  data 
with  those  of  other  workers . ^ ^ ^ ^  The  least  squares 
equations  to  all  the  data  presented  are  given  by: 
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Upper  regime 

U  =  1.573  4-  1.606  u  km/sec,  2.9  <  u  <  6.1  km/sec 
s  p  p 

S.E.(U  )  =  0.141  km/sec  for  26  points, 
s 

Lower  regime 

U  =  2.876  +  1.048  km/sec,  1.0  <  u  <  2.4  km/sec 
s  p 

S.E.(U  )  =  0.228  km/sec  for  37  points, 
s 

The  correlation  of  the  data  appears  to  be  quite  good,  particu- 

/  1  O  \ 

larly  because  the  highest  data  of  Rosenberg  '  fall  in  very 
nicely  with  the  lower  end  of  the  upper  region  described  here. 

As  in  the  dry  tuff  data,  the  problem  of  comparison  is  com¬ 
plicated  by  the  differing  initial  densities  and  the  general 
experimental  scatter.  The  point  of  initiation  of  the  high 
density  phase  appears  from  the  comparison  plots  to  occur  between 
P  *  0.25  to  0.29  Mb  and  V  -  0.290  to  0.285  cm3/gm.  The  data  do 
not  show  a  substantial  volume  decrease  but  merely  exhibit  a 
discontinuity  in  the  slope  of  the  hugoniot. 

Figure  23  displays  the  low  region  as  a  zone  bounded  by  lines 
representing  the  hugoniots  parallel  to  the  best  least  squares 
fit.  The  representation  is  given  by: 

C  2  (V  -V) 

p  _  o  o _ 

|Vo-1.048  (V-VQ)J  2 

With  Cq  from  0.268  to  0.297  cm/ysec  and  Vq  =  0.515  cm3/gm. 

The  initial  densities  for  the  various  tuffs  are  presented  in 
the  legend  in  Figure  23. 
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The  evidence  presented  here  and  for  dry  tuff  suggests  a  phase- 
transition  to  a  higher  density  form  of  some  component  of  the 
tuff.  This  component  is  most  probably  the  silicate  glass  mat¬ 
rix  of  the  tuff.  The  transition  pressure  corresponds  well  with 
the  transition  observed  by  Jones for  porous  quartz  (i.e  a 
highly  quartziferous  sandstone)  at  ^  0.260  Mb. 

Grouting  Materials  V44D  and  V6C 


The  grouting  material  tested  in  the  program  was  prepared  by 

the  Waterways  Experimental  Station  (W.E.S.),  Jackson,  Mississippi. 

Grout  V44^  is  employed  as  a  hugoniot  match  for  wet  tuff.  The 
following  is  a  brief  description  of  the  material: 


Ingredient  Wt  % 


Whitesides  epoxy  (Parts  A  and  B)  54.6 
Dolomite  sand  (passing  No.  200  screen)  19.5 
Bentonite  (aquagel)  25.9 


Details  of  preparation  and  sources  of  the  ingredients  may 
be  obtained  from  James  M.  Polatty,  Chief,  Engineering  Mechanics 
Branch,  Concrete  Division. 

3 

The  sample  received  had  a  measured  density  of  1.45  ±  0.02  gm/cm  , 

3 

slightly  lower  than  that  measured  by  W.E.S.  of  1.469  to  1.476  gm/cm 
on  a  similar  material. 

The  hugoniot  data  measured  are  listed  in  Table  VII.  Adjust¬ 
ments  to  the  shock  wave  transit  time  interval  were  made  for 
the  0.125  mm  thick  brass  rear-surface  shim.  The  time  adjust¬ 
ment  was  performed  using  the  impedance-match  method  described 
earlier  and  averaged  only  about  2%  of  the  total  time-interval. 
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Shown  in  Figures  24  and  25  are  the  data  obtained  and  the  data 

* 

for  three  tests  performed  by  W.E.S. 

A  fit  to  the  GM  data  is  given  by: 

U  =  2.998  +  1.275  u  km/sec,  1.5  <  u  <6.7  km/sec 
s  p  p 

S.E.(U  )  =  0.103  km/sec  for  7  points, 
s 

The  W.E.S.  fit  is  given  by: 

U  =  2.6b  +  1.53  u  km/sec,  1.0  <  u  <2.1  km/sec 
s  p  p 

S.E.(U  )  =  0.057  km/sec  for  3  points, 
s 

As  can  be  seen  in  Figure  26 ,  the  hugoniot  of  the  V44D  grouting 
material  lies  significantly  below  the  hugoniots  of  the  dry  and 
the  saturated  tuff.  In  the  P-u^  plane,  the  grout  is  less  stiff 
than  the  water  saturated  tuff  and  falls  below  the  dry  tuff  hug¬ 
oniot  beginning  at  'v  0.3  Mb. 

Figure  27  compares  the  hugoniot  data  for  the  W.E.S.  granite 

* 

matching  grout  V6C  and  the  GM  granite  data.  The  W.E.S.  data 
are  shown  extrapolated  to  2  Mb  and  are  seen  to  match  the  granite 
reasonably  well,  within  ru  6%  in  pressure  at  a  given  particle 
velocity . 


Private  communication,  Mr.  W.  R.  Sullivan,  Grouting  and 
Engineering  Physics  Section,  Concrete  Division,  Waterways 
Experimental  Station,  Jackson,  Mississippi  39205. 
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Schistose  Gneiss 


A  detailed  petrographic  analysis  of  the  schistose  gneiss 
tested  during  this  program  is  at  present  unavailable.  The 
following  is  a  mineralogical  description  of  the  rock  core 
used  in  the  program.* 


TABLE  VIII 

APPROXIMATE  MINERALOGICAL  COMPOSITION 
OF  SCHISTOSE  GNEISS 


Quartz 

Plagioclase 

Biotite 

Muscovite 

Clorite 

Garnet 

Magnetite  and  Pyrite 
Miscellaneous 


35%-33% 

22%-21% 

20%-15% 

10%-15% 

7% 

5% 

1% 

Trace 


"The  rock  core  is  a  homogeneous,  grey,  dense,  fresh, 
medium  grained,  quartz-biotite-plagioclase,  gneiss. 
The  mineral  assemblage  of  quartz,  plagioclase,  bio¬ 
tite,  muscovite,  and  minor  garnet  porphyroblasts , 
chlorite  and  accessory  pyrite,  magnetite,  zircon, 
apatite  and  tourmaline  as  well  as  variation  of  grain 
size  in  the  minerals  would  appear  to  indicate  that 
these  rocks  are  paragneisses ,  i.e.,  derived  by  re¬ 
gional  metamorphism  of  pre-existing  sediment." 


Petrographic  Report  No.  1M243-250,  Department  of  the  Army 
Corps  of  Engineers,  South  Atlantic  Division,  Division 
Materials  Testing  Laboratory,  Marietta  Georgia,  sampled 
by  Mr.  Canning. 
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The  hugoniot  data  are  presented  in  Table  IX.  The  initial  den- 

3 

sity  was  2.79  ±  0.02  gra/cm  .  Targets  used  were  the  rectangular 
design  described  earlier. 

The  hugoniot  data  are  shown  in  Figures  28  and  29  and  may  be 
described  by  a  least  squares  linear  equation  given  by: 

U  »  2.376  +  1.604  u  km/sec,  2.5  <  u_  <  6.1  km/sec 
s  P  P 

S.E.(Ug)  =  0,116  for  9  points. 

There  appears  to  be  a  lower  pressure  phase  below  0.3  Mb,  but 
the  two  data  points  taken  in  this  region  are  insufficient  to 
determine  an  accurate  hugoniot. 

The  data  also  exhibit  a  region  of  anomolous  behavior  in  the 
pressure  range  of  0.9  to  1.2  Mb.  The  apparent  discontinuity 
in  this  region  accounts  for  the  number  of  tests  performed  there. 
Without  further  experimentation,  the  cause  for  this  behavior 
cannot  be  determined.  If  the  data  are  analyzed  without  the  two 
points  falling  below  the  line  shown  in  the  figures  at  a 
particle  velocity  of  ^  4  km/sec  (C-1085  and  C-1086) ,  the 
linear  hugoniot  of  the  high  pressure  phase  is  given  by: 

U  =  2.363  +  1.611  u  km/sec 
s  p 

where  the  standard  error  in  U  is  reduced  to  0.098  km/sec. 

s 

Exclusion  of  these  two  shots  results  in  a  hugoniot  displaced 
0.6%  lower  in  intercept  and  0.4%  upward  in  slope. 
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Ferris  Wheel  Dolomite 

.  ( 9 ) 

The  NTS  Dolomite  is  characterized  in  Barnes.  A  detailed 

chemical  and  mineralogical  description  is  not  available,  but 

the  rock  probably  consists  of  >  95%  mineral  dolomite, 

CaMgtCO^)^,  based  on  an  x-ray  analysis  of  a  megascopically  similar 

sample  of  dolomite  from  the  Banded  Mountain  area.  The  density 

3 

of  the  samples  tested  was  2.82  ±  0.01  gm/cm  .  The  linear  hug- 
oniot  fitted  to  the  data  in  Table  X  is  given  by: 

U  =  5.196  +  1.191  u  km/sec,  2.6  <  u  <5.4  km/sec 
s  p  P 

S.E. (U  )  =  0.059  km/sec  for  4  points. 

b 

Figures  30  and  31  display  the  data  in  the  U  -u  ,  and  P-V  planes. 

®  P  (IT  Tl) 

Also  shown  are  four  data  points  obtained  by  Lombard v  '  '  for 

Area  12  dolomites  of  comparable  initial  density,  i.e.  p  ^2.84 

3  ° 

and  2.815  gm/cm  .  The  linear  hugoniot  reported  for  these  points 

is  given  by: 

U  =  5.626  +  1.016  u  km/sec 
s  p 

S.E. (U  )  =  0.048  km/sec  for  5  points. 

o 

A  linear  fit  to  the  GM  and  Lombard  data  is  yiven  by: 

Us  =  5.441  +  1.124  Up  km/sec,  0.5  <  u^  <  5.4  km/sec 
S.E.(U  )  =  0.098  km/sec  for  9  points. 


* 

Private  communications,  H.  Barnes,  USGS ,  Denver,  Colorado. 
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Banded  Mountain  Limestone 

A  petrographic  description  of  Banded  Mountain  limestone  is 
included  in  Barnes .  The  specimens  are  essentially  pure 

calcite,  96  -  98%  CaCO^  (calculated)  with  initial  densities 

of  2.702  ±  0.004  gm/cm^ .  The  line  r  least  squares  fit  to  the 
data  listed  in  Table  XI  is  given  by: 

.  U  =  4.047  +  1.391  u  km/sec,  1.0  <  u  <4.4  km/sec 
s  P  P 

S.E.(U  )  =  0.085  km/sec  for  6  points, 
s 

Figures  32  and  33  display  the  present  data  in  the  Us*-Up  and 

P-V  planes  respectively.  Included  in  Figure  32  are  the  data 

(22)  3 

of  Dremin  for  a  marble  of  initial  density  of  2.70  gm/cm 

and  the  data  of  Bass^15^  for  a  NTS  limestone  of  density 
2.66  gm/cm.  The  break  in  the  data  of  Dremin  indicate  that 
a  phase  char  occurs  in  the  marble  near  0.150  Mb;  just  below 
the  lowest  data  point.  The  lowest  point  of  Bass's  NTS  lime¬ 
stone  hugoniot  data  falls  below  a  linear  extrapolation  of  the 
higher  pressure  data  of  GM  and  Dremin  as  well  as  his  own.  This 
point  is  probably  in  a  low  density  phase  regime  similar  to  that 
reported  by  Dremin,  and  unobserved  by  the  present  work.  Evidence 
presented  here  suggests  that  only  the  top  four  data  of  Bass's 
limestone  hugoniot  are  in  the  high  pressure  phase  regime  and  we 
have  chosen  to  include  them  in  a  least  squares  analysis  for  com¬ 
parison.  The  linear  form  of  the  least  squares  analysis  of  the 
combined  data  in  the  high  pressure  regime  is  given  by: 
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U  =  4.080  +  1.386  u  ,  1.0  <  u  <  4.4  km/sec 
s  P  P 

S.E.(U  )  =  0.075  km/sec  for  10  points. 

The  combined  analysis  differs  from  the  GM  data  by  only  0.8%  in 
intercept  and  0.4%  in  slope. 


TABLE  XI 

HUGONIOT  DATA  FOR  NTS  BANDED  MOUNTAIN  LIMESTONE 


PQ  =  2.702  1 

0.004  gm/cm3 

Vq  =  0.3701 

1  0.0006  <'m3/gm 

Shot 

Us 

Shock  Velocity 
(km/sec) 

UP 

Particle  Velocity 
(km/sec) 

P 

Pressure 

(Mb) 

V 

Volume 

(cm3/gm) 

c-i2.n 

5.491* 

1.029 

0.153 

0.301 

C-1221 

6 . 820 (±0 . 013) 

1.979 

0.365 

0.263 

C-1250 

8.519* 

3.246 

0.747 

0.229 

C-1222 

9 . 436 (±0 .024) 

3.918(10.003) 

0.999 (10.003) 

0.216 (10.001) 

C-1227 

9.884* 

4.264 

1.138 

0.211 

C-1223 

10 . 318 (±0 .036) 

4.388(10.088) 

1.223(10.029) 

0.213(10.003) 

*  Single  Pin  Value 


SHOCK  VELOCITY  -  KM/SEC 
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Nevada  Desert  Alluvium 

A  chemical  and  physical  description  of  the  Nevada  desert  allu- 

(9) 

viura  is  given  m  Barnes.  The  data  are  presented  in  Table  XII 

and  Figures  34  and  35.  Probably  due  to  the  extreme  fineness 
and  homogeneity  of  the  specimens,  the  recorded  shock  velocities 
show  remarkably  little  internal  scatter  (less  than  1%  average) . 

3 

The  compaction  density  averaged  1.54  ±  0.02  gm/cm  for  these 

t  ^  \ 

tests  and  corresponded  to  the  low  density  material  of  Anderson 

(23)  3 

and  McQueen  (pQ  =  1.55  and  1.54  gm/cm  respectively). 

The  hugoniot  data  shown  in  Figure  34  does  not  appear  amenable 
to  description  by  a  single  simple  equation.  At  least  two  line 
segments  are  necessary  for  a  linear  description.  The  high  pres¬ 
sure  region  may  be  described  by: 

U  =  0.363  +  1.533  u  km/sec,  2.4  <  u  <6.1  km/sec 
»  p  P 

S.E.(Us)  =  0.132  km/sec  for  7  points. 

The  two  lowest  data  are  not  considered  adequate  in  themselves 
to  define  the  lower  hugoniot  and  no  line  or  equation  is  given. 

In  Figures  36  and  37  the  data  of  McQueen  and  Anderson  are 
compared  and  combined  with  the  present  results.  Again  the 
data  suggest  a  two  regime  hugoniot  structure.  The  high 
pressure  regime  hugoniot  is  given  by: 

Us  =  0.683  +  1.544  u  km/sec,  2.3  <  u  <  6.1  km/sec 

P  P 

S.E.(Us)  =  0.214  km/sec  for  25  points. 
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The  lower  regime  may  be  described  by: 

U  =  0.682  +  1.618  u  km/sec,  1.0  <  u  <  2.6  km/sec 
s  p  p 

S.E.(U  )  =  0.055  km/'sec  for  14  points, 
s 

The  combined  data  for  the  upper  regime  compares  with  the  GM 
data  by  raising  the  slope  by  'u  0.7%  and  the  intercept  by 
0.J  km/sec,  both  small  effects. 

f 

The  combined  data  presents  a  complicated  picture  of  the 
transition  from  the  lower  to  the  upper  regime.  Anderson's 
data  suggests  that  the  low  regime  extends  to  at  least 
0.190  Mb  and  yet  a  shot  reported  here  implies  the  minimum  of 
the  upper  regime  is  lower  than  0.145  Mb.  Whether  or  not  it 
is  possible  to  account  for  this  discrepancy  on  the  basis  of 
specimen  composition  or  experimental  technique  cannot  be  de¬ 
cided  at  present.  The  chemical  compositions  of  ours  and 
Anderson's  samples  are  nearly  identical,  but  the  mineralogy 
of  the  two  may  vary  somewhat  more . 

The  constituent  most  likely  undergoing  the  observed  transition 
between  0.145  and  0.190  Mb  is  the  quartz  component  which  amounts 
to  %  53%  in  the  initial  state.  The  large  volume  scatter  in  the 
pressure  range  of  0.20  to  0.40  Mb  makes  an  analysis  of  the  trend 
of  the  hugoniot  extremely  difficult.  With  just  the  present 
data,  the  P-V  hugoniot  appears  to  be  nearly  vertical  and  at 
higher  pressures  may  even  show  volume  expansion  due  to  shock 
heating . 


It  is  likely  that  above  0.40  Mb  the  quartz  component  is  com¬ 
pletely  converted  to  stishovite,  but  below  0.40  Mb,  it  is 
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difficult  to  determine  whether  or  not  the  quartz  component 
passes  through  a  mixed  phase  regime.  An  additional  com¬ 
plication  is  introduced  by  the  ^  36%  porosity  of  the  initial 
sample . 
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SECTION  IV 


CONCLUSIONS 


The  experimental  results  of  the  determination  of  the  hugoniot 
equations  of  state  for  several  Nevada  Test  Site  rocks  are 
summarized  in  Table  XIII.  The  present  work  has  extended  the 
ranges  of  hugoniot  equation  of  state  data  for  a  granite,  wet 
and  dry  volcanic  tuff,  dolomite,  limestone  and  desert  alluvium 
from  the  Nevada  Test  Site.  A  hugoniot.  was  determined  for  a 
schistose  gneiss  and  the  hugoniots  of  Fansteel-77  and  OFHC 
copper  were  confirmed.  In  addition,  the  range  of  hugoniot 
data  for  volcanic  tuff  impedance  matching  grouting  material 
was  extended.  In  general,  the  present  results  are  in  good 
agreement  with  data  of  other  researchers.  When  possible  a 
least  squares  fit  has  been  included  in  the  comparison  data 
for  each  reference  as  well  as  in  combination  with  the  pre¬ 
sent  results. 


TABLE  XIII  HUGONIOT  DATA  SUMMARY 
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APPENDIX 


DATA  ANALYSIS  PROCEDURES  FOR  THE 
DETERMINATION  OF  HUGONIOT  STATES 
FROM  SHOCK  WAVE  DATA 


INTRODUCTION 

The  following  is  a  detailed  description  of  the  data  analysis 
method  employed  for  the  determination  of  hugoniot  state  data 
derived  from  the  experimental  system  described  in  the  main 
body  of  this  report.  The  first  part  of  this  appendix  pre¬ 
sents  the  equations  and  assumptions  employed  in  the  analysis 
and  the  second  part  describes  a  computer  routine  written  to 
perform  the  calculations. 


Experiment  Analysis 

Figure  A-l  illustrates  the  system  schematically  with  the  follow¬ 
ing  nomenclature: 


=  impact  velocity 
=  pin  cap  thickness 
=  pin  gap 

=  target  thickness 
=  pin  position  on  line  1-A 
=  impactor  tilt 
=  shock  velocity  in  target 


cm/ y sec 


radians 

cm/ysec 
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Station 

No. 

I  1 


J  C 


t  B 


t  A 


Figure  A-l  Schematic  of  Equation  of  State  Studies  Target 


Two  cases  of  pin  correction  are  necessary;  the  first  involves 
the  "start"  pins  1  and  A  which  are  struck  by  the  impactor  directly 
The  second  case  involves  the  rear  surface  pins  B  and  C,  which 
receive  the  shock  transmitted  by  the  target  material;  and  thus 
have  closing  velocities  different  from  pins  1  and  A. 

The  times  to  pin  closure  are  given  by: 


t  =  ^1 _  +  gl 

U  .  2  u  . 

spin  pin 

t,  =  X3tan0  +  bA  +  gA 

v  U  .  2  u  . 

spin  pin 


(A-l) 


(A-2) 
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t  =  (X2  hBtanv  0>  tan  0  +  hB  ,Dm  ,  bB  .  gB 

B  - - -  _  sec(^Q)  +  ^ -  + 


2  u 

spt  pt 


(A- 3) 


= 


(X3-hctan^  8)  tan  9  + 

V  D  sec(v0^  +  U 


spt 


2  u 


pt 


(A- 4) 


where  subscripts  1,  A,  B,  C  refer  to  pin  locations,  and 
U  .  =  pin  cap  shock  velocity  induced  by  impactor 

spiii 

u  in  =  pin  cap  partical  velocity  induced  by  impactor 
ugpt  =  pin  cap  shock  velocity  induced  by  target 


u  t  =  pin  cap  particle  velocity  induced  by  target 


U 


Sp^n  and  u  -n  are  evaluated  by  the  impedance  match  solution 
for  the  direct  impact  of  the  projectile  material  into  the  pin 
material.  The  particle  velocity  derived,  from  this  analysis  is: 


B  -  V  "2 


B  -  4AC 


2A 


where  A  =  p  S  -  p  .  S  . 

op  p  opm  pin 


B=pC  +p.C.+2pCv 
op  op  opin  pin  Kop  op 


(A-5) 


(A-6) 


C  =  p  v  (C  +  S  v) 
op  op  p 


and  U  .  =  C  .  +  S  .  u  . 

spin  pm  pin  pin 
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the  impactor  hugoniot  being  given  by: 


U 

sp 


S  u 


P  P 


(A-8 ) 


For  the  pin  gap  velocity  correction  of  the  wave  interaction  be¬ 
tween  the  specimen  and  the  rear  surface  pins  on  the  impedance 
match  solution  is  employed  again.  Although  the  equations  are 
similar  to  those  already  derived,  the  hugoniot  state  of  the 
specimen  must  first  be  estimated  and  then  determined  by 
iteration. 


The  first  approximation  of  the  shock  velocity  is  found  by  the 
following  procedure:  Subtracting  Equation  (A-6)  from  Equation 
(A- 7)  we  obtain 


tan  0 


■\  [ 


‘W 


Vbl  fiLlil 

U  2  u  . 

spin  pinj 


A- 9 


Neglecting  several  small  corrections,  the  approximate  shock 
velocity  in  the  target  at  Station  B  is  given  by 


db 


'W 


V 


tan 


0 


U 

spin 


(A- 10 ) 


The  particle  velocity,  u',  at  this  shock  velocity  is  given  by 
Equation  (A- 5) ,  where: 


A  =  p  S 
op  p 


B  =  p  C  +  p  ,  D  '  4  2p  S  V 
*op  op  ot  op  p 


C  =  p  v  (C  +  S  v) 
op  op  p 


( A- 11 ) 
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The  first  approximation  of  the  shock  pressure  is  given  by: 


P'  =  p  D  '  u 1 
ot 


(A-12 ) 


The  approximate  shock  and  particle  velocities  thus  obtained 
are  used  in  the  more  exact  solution  for  shock  velocity  given 
by  Equation  (A-13) .  The  results  may  be  iterated  several  times 
until  the  desired  convergence  is  reached.  With  the  convergence 
test  set  at  0.01%  only  two  iterations  are  required. 


h  sec 

JB 


V  Q  j  fcl^ 


tan  0  jX2  -hgtan 


^-JA  + 


spin 


2  u  . 
pin 


U  ,  2 

spt 


v] 


(A-13) 


The  particle  velocity,  u  is  calculated  by  Equation  (A-5)  with: 

P 


A  =  p  S 
OP  P 


B  =  p  C  +Po.U  +  2p  Sv 
"op  op  ot  s  op  p 


C  =  p  v  (C  +  S  v) 
op  op  p 


(A-14) 


Pressure  is  calculated  from 


P  =  p  ,  U  u 
ot  s  p 


(A-15) 


where  it  is  assumed  that  P  =  0 

o 
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The  procedure  is  then  applied  to  analyze  the  data  from  pins  1, 

A  and  C,  A,  1,  and  B,  and  A,  1  and  C  using  the  appropriate  values 
of  g,  h,  b,  and  X. 

Shim  transit  time  corrections  were  made  by  adjusting  the  value 
of  b,  the  pin  cap  thickness,  to  include  the  thickness  of  the 
shim.  Since  both  shim  and  pin  cap  were  brass,  no  significant 
difference  in  shock  impedance  exists  across  the  shim-cap  inter¬ 
face. 

The  procedure  described  above  is  a  refinement  upon  the  simplest 
form  of  average  velocity  determination  (i.e.  distance/time) 
where  the  time  term  is  corrected  for  impactor  tilt.  The  above 
refinement  results  in  an  average  effect  on  Ug  of  0.2%  compared 
with  the  simplest  form.  The  overall  effect  on  a  shot  basis  has 
been  found  for  metals  such  as  copper  and  Fansteel-77  to  reduce 
data  scatter  about  a  linear  form  of  the  hugoniot  in  U  -u  . 


Data  Analysis  Program  -  SHOVEL 

The  computer  program  described  below  is  written  in  Mark  I 
Fortran  for  use  on  the  General  Electric  Timesharing  Computer 
Service.  This  Fortran  is  very  similar  to  Fortran  IV  used  on 
most  computers  and  provides  some  convenience  for  adapting  to 
teletype  I/O. 

The  program  has  two  sections  and  one  subroutine.  The  first 
portion  of  the  program  is  merely  an  operator  convenience  to 
convert  all  input  data  into  the  required  units  and  formats. 
The  latter  portion  calculates  the  shock  velocity  at  each  sta¬ 
tion.  The  subroutine  20MPUT  is  in  the  form  of  the  quadratic 

solution  that  results  in  a  meaningful  value  of  u  . 

P 
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Definition  of  Terms 


The  terms  presented  below  are  in  the  following  order: 


1.  Values  required  as  input 

2.  Intermediate  terms 

3.  Output  values 

1 .  INPUT  TERMS 

TERM 

NAME 

UNITS 

Shot  No. 

SHOT 

I 

D 

RUN 

Impactor  Material 

MATLP (2) 

H 

Target  Material 

MATLT ( 2 ) 

H 

Impactor  C 

CZERO 

km/sec 

Impactor  SU 

SP 

Impactor 

RP 

gm/cm^ 

Target  p 

RT 

gm/cnr 

Pin  1-B  ?ime 

DTB 

ysec 

Pin  1-C 

DTC 

ysec 

Pin  1  height 

HI 

in. 

Pin  A  height 

HA 

in. 

Pin  B  height 

HB 

in . 

Pin  C  height 

HC 

in. 

Pin  1  gap 

G1 

in . 

Pin  A  gap 

GA 

in. 

Pin  B  gap 

GB 

in. 

Pin  C  gap 

GC 

in . 

Pin  1-C  Spacing 

XI 

in. 

Pin  1-B  Spacing 

X2 

in. 

Pin  1-A  Spacing 

X3 

in . 

D 

H2 

D 

H3 

Polaroid  target  fiducial  spacing 

M2 

cts . 

Impactor  to  Window  #1  Ref.  Spacing 

LI 

cts . 

Impactor  to  Window  #2  Ref.  Spacing 

L2 

cts  . 

X-ray  Flash  Time  Interval 

TSI1 

ysec . 

X-ray  #1  to  Pin  C  Time  Interval 

TSI2 

ysec . 

Front  Surface  Shim  Thickness  (If  Used) 

SHIM 

in . 

Tilt  Time  Interval 

TILT 

ysec . 

Approximate  Slope  of  Target  Hugoniot 

S 

(S=1.5  if  zero) 

Pin  Combination 

COMBO 

H 

Computer  Run  No. 

RUNS 

I 

I  -  Integer,  D  -  Dummy,  H  -  Holorith  Field 
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2 .  COMPUTED  TERMS 


4 < 

TERM  NAME  UNITS 


X-ray  1  Scale  Factor 

FI 

in/ct 

X-ray  2  Scale  Factor 

F2 

in/ ct 

Window  Fiducial  Spacing 

D 

in 

Impact.  Velocity 

W 

ft/sec 

Impact  Velocity 

W1 

ft/sec 

Impact  Velocity 

V 

cm/ysec 

D 

US  TAR 

Pin  and  Po 

RPIN 

2 

gm/cm 

Pin  Material  Hugon.iot  Slope 

SPIN 

Pin  Material  CQ 

C0PIN 

km/sec 

Pin  Cap  Thickness 

B1 

in3 

Specific  Volume  (initial) 

SVOL0 

cm  /gm 

Quadratic  Solution  Term 

A 

Quadratic  Solution  Term 

B 

Quadratic  Solution  Term 

C 

Particle  Velocity  of  Pin  Cap 

UPIN 

cm/ysec 

Shock  Velocity  of  Pin  Cap 

USPIN 

cm/ysec 

Impactor  Tilt 

THETA 

radians 

First  Approximation  Us  (B  Station) 

USB1 

cm/ysec 

First  Approximation  Up 

UB1 

cm/ysec 

First  Approximation  P 

PB1 

Mb 

First  Approximation 

RT1 

gm/ciri 

First  Approximation  Sound  Velocity 

C2 

cm/ysec 

Particle  Velocity  Induced  by  Target 

UPT 

cm/ysec 

Shock  Velocity  Induced  by  Target 

USPT 

cm/ysec 

Shock  Wave  Tilt  in  Target 

THET1 

radians 

3 .  OUTPUT  TERMS 


B  Shock  Velocity 

USB 

km/sec 

B  Particle  Velocity 

UB 

km/sec 

B  Pressure 

PB 

Mb 

B  Relative  Volume 

VOL 

cm'Vgm 

Specific  Volume 

SVOL 

100 


-  ..  -,vi  -<rv 
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The  term  labeling  is  essentially  identical  for  the  "C"  station 
calculations  with  the  output  being: 


C  Shock  Velocity 

use 

km/ sec 

C  Particle  Velocity 

UC 

km/sec 

C  Pressure 

PC 

Mb 

C  Relative  Volume 

V0L1 

3 

C  Specific  Volume 

SV0L1 

cm  /gm 

The  "D"  station  calculation  employs 

the  C  station  approximate 

velocity  and  corrects  for  the  time 

interval  difference 

between 

B  and  C.  The  output  terms  are: 

Shock  Velocity 

US 

km/sec 

Particle  Velocity 

U 

km/sec 

Pressure 

P 

Mb 

Relative  Volume' 

VOL  2 

Specific  Volume 

SV0L2 

2 

cm  /gm 

Additional  output  terms: 


D  Station  Calculation  of  Shock  Tilt  THET3  Radians 

Calculated  Impactor  Tilt  Time  Interval  TILT1  ysec 

%  Difference  of  (B-C) /C  Shock  Velocities  DIFF  % 


The  Mark  II  Time  Sharing  System  allows  chaining  of  programs  and 
SHOVEL  is  continued  as  SHOVEN  and  an  external  file  INFORM  is  used 
as  the  data  source.  Up  to  seven  files  may  be  linked  and  called 
in  process. 
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SHOVEL 

100  COMMON  U,A#B*C 

1  10  REAL  DIMENSION  MATLPC2).  MATLT(2> 

120  INTEGER  SHOT*  RUN# COMBO#  RUNS 
1  30C 

140  200  READC  1  #  40)  SHOT#  RUN#  MATL PC  1  )  *  MATL  PC  2  ) »  MATL  TC  1  )  »  MATLTC  2  ) 

1 50  READC 1 )CZERO#  SP*  RP*  RT#  DTB# DTC# H 1  *  HA* HB# HC* G I »  GA.GB*  GC*X1*  X2*  X3* 

1 60  ♦H2#H3#M2#L1#L2*  TSI 1*  TSI2* SHIM#  TILT#  S 
1  70C 

180  40  FORMATC  5X#  15#  12#  4A< 

1  90C 

200  PRINT# "ENTER  PIN  COMBINATION  <1  START  =  1#  A  START  *  A)  AND  RUN  NO." 

2  10  INPUT  30*  COMBO#  RUN  S 
220C 

2  30  30  FORMAT<2Al> 

2  40C 

250CMPTV  COMPUTE  TILT  AND  VELOCITY 

252  PRINT  61#CZER0»SP»RT*RP*DT8*DTC#HB#HC#G1#GA»GB,GC#X1#X2#X3# 

254  +TILT#  S 

256  61  FORMAT<////6F5.3#2F7.5#4F6.4//5F6.4///> 

260C 

270  IFCL1 )46*48#46 

280  46  Fl=. 4756E-03lF2=.4771E-03lD=l 1.9951 

290  IFCL2141#  42#  41 

300  41  D 1  * (L 1  *F1 +L2+F2+D)/  1  2*  0 

310  VEL1=*CD1/TSI  1  )*1.0E  +  06 

320  42  D2*CL1*F1+D*M2*F2*.9537*SHIM>/12.0 

3  30  VEL2*  C  02/  <  TSI  2-  DTC  )  )  *  1  •  0E+06 
340  IF(L2)43»  44#  43 

3  50  44  GO  TO  45 

360  43  W=  VEL 1 1  W 1  *  VEL2 

3  70  GO  TO  49 

380  45  W*VEL2l  W1*0.  0 

3  90  GO  TO  49 

400  48  READC 1)W 

410  W1=0.0 

420C 

430CTC  CONVERT  INPUT  TO  CGS 

4  40C 

450  49  V=  C W/3»  281 >* 1 • 0E+02 

460  CZERO*CZERO* 1 • 0E+05 1 USTAR*  USTAR* 1 »  0E  +  05 

4  70  DTB* DTB* 1 • 0E- 06 f  DTC* DTC* 1 • 0E-061  TIL T= TILT* 1 . 0E- 06 
480  Hl*Hl*2»54l  HA* HA* 2. 541HB=HB*2. 541 1HC=HC*2.54 

490  RPIN  s8.41lSPIN=1.42lC0PIN=3.8E  +  05lBI*5.0SE-03 

5  00  G1  =  G1 *2. 541 GA*GA*2. 541GB* GB*2 . 541 GC=GC*2.  54 
510  X1=X1*2.54!X2=X2*2. 54lX3=X3*2. 54 

520  THET3=0.  01  SVOL0*  1 . 0/RTl  SVOL=0.  01  SVOL  1*0.  0 
5  30C 

5  40CSV  COMPETE  SHOCK  VELOCITY 
5  50C 

560  IF  CDTBll »* 1 1* 1000 
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SHOVEL  CONTINUED 


5  70  II  UB-0.0lPB«0*0lUS8=0.0lVOL»0*0 
5  80  GO  TO  5 

5  90C 

600  1000  B«RP*CZERO+RPIN*C0PIN+2*RP*SP*V 

6  10  A-RP*SP-RPIN*SPJNlC»RP*V*CCZEROt-SP*V> 

620  CALL  COMPUT 

630  UPIN-U1USPIN-C0PIN+SPIN*UPIN 

6  40  THETA- C  V/X3)  *(  TIL  T-CGA-Gl  >/<2*UPIN>  > 

650  USBI* (H8-H1 >/<DTB-X2*THETA/ V> 

660  1001  8- R P*C ZERO +RT* USB  1 ♦2*RP*SP*Vl A-RP*SP 
670  C*  RP*V*  (C  ZERO*  SP*V) 

680  CALL  COMPUT 

690  U81-UlPBl-RT*USBl*UBl 

700  RT1«CRT*USBI )/CUS81-UBl ) 

7 10  C2-USB1 *SQRT( • 49+  <  <  USB 1 • UB1 >/USBl >**2> 

720  I  F<  S)  50*  50#  5  1 

7  30  50  S-  1  •  5 

7  40  51  B-RPIN*C0PIN+RT1*C2+2*RT1*S*UBI 1A-S*RT!-SPIN*RPIN 
7S0  C*  S*RTI  *U81  **2*RT1  *C2*UB1  +P8 1 
760  CALL  COMPUT 

7  70  UPT-Ul  USPT-C0PIN  +SPIN4UPT 
780  THETl-<  USB 1/V>* THETA 

790  DELT1-OTB-  THETA*<  X2-  CHB-H 1  >  *THET1  )/V 

8  00  DELT1*DELT1+B1/USPIN+G1/<2*UPIN>-8I/USPT-GB/(2»UPT> 

8  10  USB-CHB-MI )*SQRT(  THET! **2+1 )/DELTl 

820  3  I  F(  *  0001 - ( ABSC  C USB- US81 >/USB>  > ) I  002# 1 003# 1 003 
830  1002  USB 1- USB 
8  40  GO  TO  1001 

850  1003  BaRP*CZERO+RT*USB+2*RP*SP*VlA*RP*SPlC-RP*V*CCZERO*SP*V> 

860  CALL  COMPUT 

8  70  UB- Ul PB- RT*USB* UBl VOL- 1 - U8/ USB 
8  80  SVOL0*1.0/RTiSVOL*VOL/RT 

8  90C 

900  5  IF(DTC>6*  6»  7 

9  10  6  USC*  0.  01  UC*  0.  01  PC*  0.  01  VOL  1*0«  01  TILT*  TILT*!*  0E+061  US=0*  0 1 U*  0.  0 
9  20  TI L  T1 * 0*  0 1 USB*  USB* 1  *  0E- 05 1 UB*  UB* 1  *  0E- 051 PB*  PB*  I  *  0E- 1 2 

930  P-0.01  VOL2-0*  0 

940  SVOL  I  *0.  01  SVOL2*  0.01DIFF-0.  0 

9  50  GO  TO  2 

9  60C 

970  7  B»RP*CZERO+RPIN*C0PIN+2*RP*SP*V 
9  80  A«RP*SP-RPIN*SPINlC»RP*V*<CZERO  +  SP*V> 

9  90  CALL  COMPUT 

1  000  UPIN-U1USPIN*C0PIN+SPIN*UPIN 
1  010  THETA- <V/X3)*< TILT- < GA- G 1 )/<2*UPIN) ) 

1 020  USC1*<HC-H1 >/(DTC-Xl * THETA/ V) 

1030  1004  B*RP*CZER0+RT*USCl+2*RP*SP*Vl A-RP*SP 

1040  C«RP*V*CCZERO  +  SP*VJ 

1050  CALL  COMPUT 

1060  UCl-Ul  PC1*RT*USC1*UC1 
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MSL-68-15 


SHOVEL  CONTINUED 


1  070  RT2=<RT*USC1 J/CUSC1-UC1 > 

1  080  C3*U$C1*SQRT(. 49  +  C  <USC1-UC1 J/USC1 >**2> 

I  090  B*RPIN*C0PIN+RT2*C3+2*RT2*S*UC1 * A* S*RT2- SPIN*RPIN 
1  100  C»S*RT2*UC1**2+RT2*C3*UC1+PCI 
1110  CALL  COMPUT 

1 i20  UPTl*Ul  USPT1=C0PIN+SPIN*UPT1 

1  130  THET2*<USC1/V)*THETA 

1  140  DEL T2* DTC- THETA* ( X I - < HC- H ! )*THET2)/V 

1  150  DELT2=«DELT2+B1/USPIN*G1/<2*UPIN>-B1/USPT1-GC/C2*UPT1  > 

1 160  USC-CHC-Hl )*S0RT(THET2**2+I )/DELT2 

1170  4  IFC.0001-C ABS<<USC-USC1 >/USC> > > 1 005# 1 006* 1 006 

1180  1005  USCI*USC 

1190  GO  TO  1004 

1 200  1006  8«RP*CZER0+RT*USC+2*RP*SP*VIA»RP*SPIC=RP*V*CCZER0+SP*V) 
1210  CALL  COMPUT 

1220  UC»UlPC«RT*USC*UC*  VOL  I* 1-UC/USC 
1  230  UB*UB* 1  •  0E- 5l US8=US8* 1  •  0E- 05  >  PB*  PB*  1  • 0E- 1 2 
1  240  UC=  UC* 1 . 0E- 05 1 USC* USC* 1 • 0E- 05 *  PC* PC* 1 . 0E- 1 2 
1250  DELT*  (  OELT1  -  DEL  T2)  +  l*0e>06 
1260  TILT=TILT*1. 0E+06 

1  270  TILTI*(DT8-DTC)*X3/<X2-XI >*! *  0E+06 
1  280  DIFF*0.0ISVOL1*VOL1/RT 
1  290C 
1  300C 

1310  I  F<  DTB  )  7 0*  7 0#  2 1 

1320  70  TILT1*0.0IU=0.0IUS*0*0IP*0.  0IVOL2*0.  0 
1  330  GO  TO  2 
1  340  21  US1=(HC-H1  >/DTC 

1 350  B=RP*CZER0*RT*USl+2*RP*SP*VlA=RP*SP 
1360  C«RP*V*<CZERO*SP*V> 

1  370  CALL  COMPUT 

1  380  Ul*UlPl  =  RT*USl*Ul 

1390  1007  RT3*<RT*US1 >/<USl-Ul > 

1  400  C4=>US1*S9RTC  .  49*C<US1-UI  )/USl >**2> 

1410  I FC  S)22#  22*23 
1  420  22  S= 1 . 5 

1430  23  9=RPIN+C0PIN+RT3*C4+2*RT3*S*U1 1 A=S*RT3-SPIN*RPIN 
1440  C=S*RT3*Ut**2+RT3*C4*Ul+Pl 
1 450  CALL  COMPUT 
1  460  UPT*Ul  USPT«*C0PIN  +  SPIN*UPT 

1 470  DTILT*DTB- (HB-H1 >/US 1 +B l /USPIN+G I / C 2*UPIN > -B 1 /U3PT- GB/ ( 2*UPT) 

1  430  THET3=V*DTILT/X2 

1 490  SHTIM*DTC-THET3*(X1-(HC-H1 > *US 1  * THE T3/ V) / V+B 1 /USPIN+G 1 / ( 2* UPI N > 
1500  +-B1/USPT-GC/(2.*UPT) 

1510  US3  <HC~H1 >*SQRT< 1. 0+<USI*THET3/V>**2>/SHTIM 
1  520  B*RP*CZER0+RT*US+2*RP*SP*VJ A*RP*SP 
1530  C=RP*V*<CZERO+SP*V> 

1540  CALL  COMPUT 
1  550  U=*UJ  P=RT*US*U 

1 560  I F( #  0001 - ( ABS<  <  US- US  1 )/US> ))  1003*  1009*  1009 
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MSL-68-15 


SHOVEL  CONTINUED 


1  570  1008  US1  =  USIU1  =  U»P1  =  P 
1530  GO  TO  1007 

1  590  1009  V0L2=  1  -  U/USl  U*  U*  1  »  0E-  05  J  US*  US*  1  •  0E-05*  P*  P*  1  »  0E-  1  2 
1  600  SVOL  I  *  VOL  1  /RTl  SV0L2= V0L2/RT 
1610  DI FF= ( ABSC  C  USB- US) /USC ) ) *  1 00 • 0 
1 620  SUSE  SHOVEN 


SHOVEN 


1700C  OUTPUT  STATEMENTS  FOR  SHOVEL 
1  710C 

1720  2  PRINT  60*  SHO  T*  RUNS*  COMBO,  MATL  TC  I  >»MATLTC2>*  MATL  PC  l  )  *  MATL  PC 2) * 

1  733  +RT*  RP*  SVOL0*  W*  THET3*  V*TILT*W1*TILT1*DI  FF* 

1 740  +038*  U8* PB» VOL*  SVOL*  USC*  UC*  PC*  VOL  1  *  SVOL 1  *  US*  U, P, V0L2*  SVOL 2 
I  7S0C 

1760C  OUTPUT  FORMAT 
1  7  70C 

1780  60  FORMAT*//////// 10X*  9HSH0T  NO.  *  I  5,  1  0X*  8HRUN  NO.  ,A!// 

1  790  +  1  0X*  1  7HPIN  COMBINATION:  *A1*IX*7H-  START,// 

1800  +5X*  14HTARGET  MATL.-  *  2A6* 

1810  +4X* 13HPR0J.  MATL.-  * 2A6//8X, 9HDENSI TY-  *F6»3*  IX*  5HGM/CC*  8X, 

1820  ♦9HDENSI  TY-  *  F6.  3*  1  X,  5HGM/CC//8X*  8HV0LUME-  *  F6*  4*  1  X*  5HCC/GM* 

1830  +/////5X* 18HPR0JECTILE  VEL.  *  *EI2.5* 

1840  ♦1X*6HFT/SEC,4X*  13HSH0CK  TILT  »  *  F9. 6//23X*  E 1 2*  5*  IX*  6HCM/SEC* 

1  850  +4X*  1  3HTILT  MEAS.  =  *  F9.  6//23X*  E  1  2«  5,  IX*  6HFT/SEC*  4X* 

1860  -M3HTILT  CALC.  *  *  F9.  6///5X*  2 1 HSHOCK  VEL.  SCATTER  =  *  F7.  3*  1  X*  1 H#* 
1 870  +//////2X* 4HSTA. 4X* 1 0HSHOCX  VEL.  4X» 1 0HPART*  VEL-  7X*  8HPRESSURE 
1  880  ♦  4X»  4HV/V05X*6HVOLUME//10X*6HKM/SEC,8X,  6HKM/SEC*  1  1X*8HMEG^  'ORS* 
1890  ♦ 1 3X*  5HCC/QM//// 

1900  +2X*  1  KB*  6X*  F8*  4*  6X*  F8®  4*8X*F10.6*  1X*F9.5*  1X*F8.5// 

1910  +2X*  IMG*  6X* F8«  4*  6X*F8.  4*  8X*F10.  6*  I  X.  F9 . 5*  1  X* F8. 5// 

1920  +2X*  1HD*6X*F8.4,6X,F8.4*8X*F10.6*  IX.F9.5*  1  X,  F8.  5////// ) 

1 930  GO  TO  200 
1  940  END 

1  950  SUBROUTINE  COMPUT 

1960  COMMON  U*A,8*C 

1  970  U=*  <  B-  SORT*  B**2-  4*A*C))/*2*A> 

1980  RETURN 
1 990  END 

2000  SFILE  INFORM 
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